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Chapter 



Preface 



Time, and equivalently frequency, is the most precisely measured physical 
quantity. It is therefore inevitable that virtually all domains of engineering 
and physics need reference oscillators. The oscillator noise can be decom- 
posed into amplitude noise and phase noise. The latter, far more important, 
affects timing, for it is related to precision and accuracy of measurements. 

The oscillator, inherently, turns the phase noise of the internal parts 
into frequency noise. This is a necessary consequence of the Barkhausen 
condition, which states that the loop gain must be of one, with zero phase, 
for stationary oscillation. There follows that oscillator phase noise, which 
is the integral of the frequency noise, diverges in the long run. This phe- 
nomenon is often referred to as the "Leeson model" after a short article 
published in 1966 by David B. Leeson |Lee66j . and called Leeson effect 
here, in order to emphasize that it is far more general than model. In 2001, 
David B. Leeson received the W. G. Cady award of the IEEE International 
Frequency Control Symposium "For clear physical insight and model of the 
effects of noise on oscillators" . 

Since spring 2004, I had the opportunity to give some seminars on noise 
in oscillators at the Jet Propulsion Laboratory, at the IEEE Frequency 
Control Symposium, at the FEMTO-ST Laboratory, and at the Universite 
Henri Poincare. These seminars had the purpose to provide a tutorial, as 
opposed to a report on advanced science, addressed to a variety of people 
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including technicians, PhD students, and senior scientists. This monograph 
derives from these seminars, and from numerous discussions with coheagues. 

The topics covered can be divided into three parts. Chapter 1 addresses 
language and general physical mechanisms. Chapter 2 aims at understand- 
ing the inside of commercial oscillators through the analysis of the speci- 
fications. Chapter 3 and 4 focus on the use of the Laplace transform to 
describe the oscillator and its phase noise. 

Nancy, Feb 23, 2005 



Enrico Rubiola 



Chapter 1 



Heuristic approach to the Leeson effect 



1.1 Phase noise fundamentals 

This introductory section provides a summary about phase noise and of 
its properties. The material is available in many classical references, such 

|H,ut78i RTTmnl IvM |vii99 . 

The quasi-perfect sinusoidal signal of oscillators is modeled as 

v{t) = Vo[l + a{t)] cos[27rzyot + <f{t)] , (1.1) 

where fo is the carrier frequency; the random variables a{t) and ip{t) are the 
fractional amplitude noise and the phase noise, respectively. The physical 
dimension of ip{t) is rad, a{t) is dimensionless. 

It is sometimes convenient to rewrite the signal p.l|) in the equivalent 
Cartesian form 

v{t) = Vq cos(27rfoi) + Vc{t) cos(27rfot) — Vs{t) cos(27r;/ot) . (1-2) 
In low noise conditions (|a| <^ 1 and \(p\ <^ 1), it holds that 

a(t) = ^ and ^(0 = ^- (1-3) 

In the absence of noise, the spectrum of v{t) is a Dirac Vq — uq) 
function. Noise broadens the spectrum. Most of the art of measuring the 
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oscillator noise is related to the ability to measure extremely narrow-band 
signals, for the radiofrequency spectrum turns out to be a poor tool. The 
oscillator noise is better described in terms of the power spectrum density 
S{f) of the amplitude and phase noise, thus Sa{f) and S^{f), as a function 
of the Fourier frequency /. Only phase noise is analyzed here. Nonetheless, 
one should be aware that the effect of amplitude noise may not be negligible, 
and that the resonant frequency of some resonators may be affected by the 
amplitude. 

The physical unit of S^{f) is rad^/Hz. Phase noise spectra are (almost) 
always plotted on a log-log scale. The technical unit "decibel", 5dB = 
10 log^oC*?), is commonly used. Manufacturers prefer the quantity C{f) 
(pronounce 'script-ell') to S^{f). In physics and mathematics S^{f) is 
preferred. £(/) and S^{f) are equivalent since C{f) is now^ defined as 



£(/) is always given in dBc/Hz, which stands for dB below the carrier in 
1-Hz bandwidth. In decibels, £(/) = 5'^(/) — 3dB. 

A model that has been found useful in describing the oscillator noise 
spectra is the power-law 



i=0 

Table ITTTI shows the phase noise terms of (|1.5j) . If needed, the sum (jl.sp 

may be extended adding additional negative terms. 

Two other quantities are often used to characterize the oscillator noise, 



(1.4) 



-4 




(1.5) 



x{t) 



phase time 



(1.6) 



yit) 



2'kvq 



fractional frequency fluctuation . 



(1.7) 



^Formerly, £(/) was defined as the single-sideband noise power in 1 Hz bandwidth 
divided by the carrier power. This definition has been superseded by 11.41 because it was 
ambiguous when amphtude noise and phase noise have not the same spectrum. 
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Table 1.1: Noise types, power spectral densities, and Allan variance. 



noise type 




SyU) 


Sip Sy 




mod(T^(r) 


white if 


bo 




h2 = ^ 
^1 = ^ 

ho = '^ 


OC T"^ 


OC 


flicker if 




hif 


OC T"^ 


OC 


white / 




ho 




\hor~^ 


flicker / 


b-sf-' 




h-1 = ^ 

h 2 ~ 


21n(2) /i_i 


i ln(2) 


rand, walk / 


b-4f-^ 


h^2r' 


— h-2T 


— h-2T 



The phase time (fluctuation) x{t) is the phase fluctuation ip{t) converted 
into time, and measured in seconds. The fractional frequency fluctuation 
y{t) is the instantaneous frequency fluctuation normalized to the nominal 
carrier frequency uq. y{t) is dimensionless. The power spectral densities 
are 

S.U) = \s^{f) (1-8) 

^0 

Sy{f) = ^SM) ■ (1-9) 

Sy{f) [Ea. (|1.9jl ] is obtained from the definition p.7jl using the property 
that the time-domain derivative maps into a multiplication by juj = j27r/ 
in the Fourier transform domain, thus by (j^ = 47r^/^ in the spectrum. 

The power-law model applies to Sx{f) and Sy{f). The coefficients of 
Sy{f) are denoted by hi in the literature, hence 

-2 

Sy{f)=Y.hif . (1.10) 

i=2 



Table ITTTI helps in conversions between Sip{f) and Sy{f). 
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Figure 1.1: Power-law, spectra, and Allan variance. 



Another tool often used in the oscillator characterization is the Allan 
variance cry(r), as a function of the measurement time r. The Allan vari- 
ance is always estimated by averaging. Given a stream of M data y, each 
representing a measure of the quantity y{t) averaged over a duration r 
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arg(Ap) 
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Figure 1.2: Basic feedback oscillator. 



ending at the time t = kr, the estimated Allan variance is 



1 



2(M - 1; 



A/-1 

E 

fc=i 



Vk+l ~ Vk) 



(1.11) 



Table ITTTI provides some conversion formulae to calculate o'y(r) from Sy{f). 
It is important to understand that c^Ct) can always be calculated from 
Sy{f), but the inverse is not free from errors |Gre98j in the general case. 
The modified Allan variance, not analyzed here, is also commonly used. 
Figure ITTTI provides a summary of the power-law spectra and Allan variance. 



1.2 Oscillator fundamentals 

The basic feedback oscillator (Fig. \1.'2\\ is a loop in which the gain A of the 
sustaining amplifier compensates for the loss [gain f3{uj) in the figure] of the 
resonator at a given angular frequency frequency ujq. The condition for the 
oscillation to be stationary, known as the Barkhausen condition, is 

A(3{uj) = 1 Barkhausen (1.12) 

\Af3{u;)\ = 1 (1.13) 
aig[A(3{uj)] = (1.14) 



at u) = ajQ. 
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amplifier 



resonator 



amplifier 



resonator 



R<0 





I 



V 



Figure 1.3: Negative-resistance oscillator. 



The unused input (0 V) in Fig. 11.21 serves to set the initial condition 
that triggers the oscillation, and to introduce noise in the loop. 

It is often convenient use a constant-gain amplifier (A is independent of 
frequency), and a bandpass filter as /3 = /3(a;) in the feedback path. Some 
small frequency dependence of the amplifier gain, which is always present 
in real- world amplifier, can be moved from A to P = (5{io). The function 
(3{uj), still unspecified, is described graphically in Fig. 11.21 

The model of Fig. 11.21 is general. It applies to a variety of systems, 
electrical, mechanical, and others. A little effort may be necessary to iden- 
tify A and /5. If for example the resonator is a two-port microwave cavity 
connected to an amplifier in closed loop, matching it to Fig. 11.21 is trivial. 
A less trivial example is the negative-resistance oscillator shown in Fig. 11.31 
The feedback function P{u}) is the resonator impedance Z{u;) = thus 
I{uj) is the input and V{u!) the output. The resonator impedance is a com- 
plex function of frequency that takes a real value (a resistance) at lo = ujq. 
A negative conductance G plays the role of the amplifier. We match the 
oscillator of Fig. 11.31 to the general scheme fFig. 11.21) observing that the 
controlled current generator is a transimpedance amplifier that senses the 
voltage V across the resonator and delivers a current / = GmV. The game 
of signs deserves some attention. The condition \AP\ = 1 requires that 
Gm > 0. The sign of the current can follow two conventions, in a generator 
the current is positive when it exits, in a load the the current is positive 
when it enters. Interpreting the controlled generator as a resistor, the sign 
of the current is to be changed. Thus Gm = —G. 

Oscillation starts from noise or from the switch-on transient. In the 
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spectrum of such random signal, only a small energy is initially contained 
at loq. For the oscillation to grow up to a desired amplitude, it is necessary 
that |^/3(u;)| > 1 at uj = ujq for small signals (Fig. 11.41) . such condition, 
oscillation at the frequency ujq that derives from arg[^/3(u;)] = rises ex- 
ponentially. As the oscillation amplitude approaches the desired value, an 
amplitude control (not shown Fig. II. 4() reduces the loop gain, so that the 
loop reaches the stationary condition A(3(ij) = 1. The amplitude can be 
stabilized by an external AGC (automatic gain control), or by the large- 
signal saturation of the amplifier. Figure [131 shows the effect of saturation. 
When the input amplitude exceeds the saturation level, the output signal is 
clipped. Further increasing the input level, the gain decreases at the funda- 
mental frequency wq, and the excess power is squeezed into the harmonics 
at frequencies multiple of loq. 

In summary, it is important to understand that in real-world oscillators 

1. it is necessary that |^/?(w)| > 1 for small signals, 

2. the condition |j4/3(w)| = 1 results from large-signal gain saturation, 

3. the oscillation frequency is determined only by the phase condition 
arg[Al3{uj)] = 0. 
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v(t) clipping S(co) excess power goes 

k A to higher frequencies 




Figure 1.5: Time-domain clipping results in power squeezed to higher har- 
monics. In the example shown clipping is almost symmetric, for the odd 
harmonics are privileged, but the even harmonics are still present. 



If a static phase ip is inserted in the loop (Fig. II. HI right V the Barkhausen 
phase condition becomes arg/?(a;) + ^ = 0. Hence the loop oscillates at the 
frequency 

0^0 + AliJ at which arg/3(a;) = —ip . (1-15) 

The effect of on the oscillation frequency is obtained by inverting Eg. 11.151 
Within the accuracy of linearization, it holds that 

AiO = . (1.16) 

f;arg/3(..) 

If the resonator is a simple circuit governed by a second-order differential 
equation with low damping factor (i.e., large merit factor Q), in the vicinity 
of the resonant frequency it holds that ^ aig P{uj) = — Thus 

Auj Au ijj ^ Auj I , , 

= = — — resonator, for <^ —— . (l-l') 

ujQ vq 2Q loq 2Q 

For reference, ip > means that the loop leads in the time domain, conse- 
quently the oscillator is pulled to a frequency higher than the exact reso- 
nance. 



1.2.1 Pulling the oscillator frequency 

There exist (at least) two ways to tune an oscillator to the desired frequency. 
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static phase 





arg(AP) 




CO 



^ tuning range 
arg (AP)+\tf=0 



oscillation cOg+Aco 



Figure 1.6: Tuning the oscillation frequency by insertion of a static phase. 



Introduce a static phase shift in the loop The first method consists 
of introducing a static phase lead or lag "0 in the loop, as in Fig. 11.61 
Oscillation is ruled by the Barkhausen condition (|1.12j) . with saturated 
amplitude. Therefore, the oscillator tuning range the frequency range in 
which 



> 1 (smah signal) 



(1.18) 



In this region, the gain can be reduced by saturation and the phase de- 
termines the oscillation frequency. Out of this range, the response to a 
perturbation decays exponentially, hence no oscillation is possible. This 
method is used in microwave oscillators, where a phase shifter is used to 
set the static phase ip. 



Change the natural frequency of the resonator The second method 
consists of pulling the natural frequency of the resonator by modifying the 
parameters of the resonator differential equation. The adjustment circuit 
is no longer distinct from the resonator, and there is no reason to introduce 
the static phase tp. This method is often used in quartz oscillators, where a 
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quartz resonator 
equivalent 



electrode 



capacitance 

-AV 



external 
frequency adjustment 



Figure 1.7: Typical tuning scheme for quartz oscillators. 



variable capacitor is used to alter on the resonator natural frequency (Fig. 



The effect on phase noise Analyzing noise, the two methods are quite 
different. The static phase tp, inherently, increases the noise bandwidth of 
the resonator. Conversely, the reactance used to pull the natural frequency 
of the resonator has not such intrinsic effect. Nonetheless, in practice the 
resonator noise bandwidth still increases because the additional loss intro- 
duced by the external reactance reduces the merit factor. More details are 
given in Section [3.4. II 



1.3 The Leeson formula 

Let us consider an oscillator in which the feedback circuit P is an ideal 
resonator, free from frequency fluctuations, with a large^ merit factor Q. 
The resonator relaxation time is 

T = —To = = — (1.19) 

Let us then replace the static phase of Eq. (|1.16l) with a random phase 
fluctuation 'ip{t) (Fig. I1.H|1 that accounts for all the phase noise sources in 

^Strictly, only Q S> 1 is necessary. 
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random phase 



out 



Figure 1.8: The phase noise of the amphfier and of all other components of 
the loop is modeled as a random phase "0 at the input of the amplifier. 



the loop. There follows an oscillator output signal of the form 

v{t) = Vo cos[u;ot + ^(t)] , (1.20) 

where ip{t) is the effect of ip{t). We analyze the mechanism with which the 
power spectrum density of ip is transferred into if. 

For the slow components of ^{t), slower than the inverse of the relax- 
ation time, ip can be treated as quasi-static perturbation. Hence 

A^=^V'(t) (1.21) 

and 

5a.(/)= (^)'5^(/) . (1.22) 

The instantaneous output phase is 

ip{t) = 2vr y (Ai/) dt . (1.23) 

The time-domain integration maps into a multiplication by i in the Fourier 
transform domain, thus into a multiplication by ™ spectrum. 

Consequently, the oscillator spectrum density is 



sAf) = 72 (^^ ) sAf) ■ (1-24) 
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For the fast fluctuations of ip, faster than the inverse of the relaxation 
time, the resonator is flywheel that steers the signal. Loosely speaking, 
it is open circuit for the phase fluctuation. The fluctuation ■ip{t) crosses 
the amplifier and shows up at the output, without being fed back at the 
amplifier input. No noise regeneration takes place in this conditions, thus 
ip{t) = tp{t), and 



(1.25) 



Under the assumption that there is no correlation between fast and slow 
fiuctuations, we can add the effects stated by Equations (|1.24|) and H1.25|) 



P {2Qj 



Leeson 
formula 



The above can be rewritten as 



1 + 



fl 
P 



SM) > 



(1.26) 



(1.27) 



where 



f = ^ 
2Q 



2ttt 



Leeson frequency 



(1.28) 



is the Leeson frequency. By inspection on Eq. ()1.27|) . the oscillator behavior 
is that of a first-order filter with a perfect integrator (a pole in the origin in 
the Laplace transform domain) and a cutoff frequency /x, (a zero on the real 
left-axis), as shown in Fig. 11.91 The filter time constant is the relaxation 
time T of the resonator. 

It is to be made clear that Eq. (|1.26jl . and equivalently Eq. (|1.27j) . 
accounts only for the phase-to- frequency conversion inherent in the loop. 
The resonator noise is still to be added for the noise spectrum to be correct. 

To summarize, the Leeson effect [Eq. (|1.26j) - (|1.27j) ] consists of the mul- 
tiplication by of the amplifier phase noise spectrum below the Leeson 
frequency /i = |q- This behavior is quite general since the amplifier noise 
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Ml 



log-log 
scale 



s„,(f) 



1 + 



S^(f) L fi\2Qj 



Leeson effect 



Figure 1.9: The Leeson-efFect transfer function. 




oscillator noise 
amplifier noise 



Ojo \ -2 white frequency 



white phase bg 



amplifier white 



Figure 1.10: The Leeson effect, shown for a typical microwave oscillator. 



is still unspecified. Figure II.IUI shows the Leeson effect in a typical case 
(microwave oscillator) , in which the amplifier shows white and flicker phase 
noise. 

The formula (|1.26|) was originally proposed by David B. Leeson |Lee66j 
as a model for short-term frequency fluctuations, inspired to the magnetron 
for radar applications. It was perfectly sound from this standpoint to con- 
sider the cavity ideally stable in the short term (1 ms corresponds to a 
round-trip of 300 km), and to ascribe all the noise to the amplifier, which 
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relies on an electron beam. 



1.3.1 Delay- line oscillator 

The frequency reference can be a delay line instead of the resonator, as 
shown in Fig. 11.111 In the frequency domain, the delay line is described 
by P{lo) = e~^^'^ . Thus the loop can sustain any oscillation frequency for 
which arg/3(w) = 0. A selector circuit, not shown in Fig. II. ill is therefore 
necessary to select a specific oscillation frequency ojq. 

The Leeson effect is derived in quasistatic conditions from Eq. H1.16() 

Aw = . (1.29) 

^arg/3H 

In the case of the delay line it holds that ^ arg/3(u;) = — r. Interestingly, 
the delay line is equivalent to a resonator of resonant frequency vq and 
merit factor 

Q = TTl^QT , (1.30) 

thus 

/l = — (1.31) 

vrr 

to the extent that it has the same slope ^ arg/3(u;) of the resonator. 
For slow fluctuations it holds that 

Ai/ = ^ / « /l , (1.32) 

and therefore 

SAf) = ^2 ^ S^(f) f«fL- (1-33) 

The noise propagation oifast phase fluctuations (/ > fi) from the amplifier 
input to the oscillator output is far more complex. Chapter |1] is devoted to 
this topic. 
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Figure 1.11: Basic delay- line oscillator. 



1.4 Amplifier noise 
1.4.1 Additive white noise 

The amplifier noise is described in terms of the noise temperature Ta defined 
as follows. When the amplifier is input-terminated to a resistor at the 
temperature Tq, the equivalent spectrum density at the amplifier input is 
= k{Ta + To). The amplifier noise is therefore 

Na = kTa . (1.34) 

The spectrum density of the equivalent input noise can also be written as 

Na = FkTo , (1.35) 

The above Eq. p.35|) defines the noise figure F. By equating ()1.34|) and 
Xrm . we find 

Na = {F-l)kTo, (1.36) 



and 

F = 



Tg+Tp 

To 



(1.37) 
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The unambiguous definition of F requires that the temperature Tq is spec- 
ified. The standard value is 290 K (17 °C). Accordingly, it holds that 
kTo = 4x10"^^ J, that is, -174 dBm in 1 Hz bandwidth. 

When amplifiers are cascaded, the input noise — 1 of the n-th am- 
plifier is divided by the power gain nm^\ of the preceding amplifiers. 
The noise figure of the chain is given by the Friis formula |Fri44j 

^ = ^^ + %r + 'lp|' + --- formula (^'^^^ 

The typical noise figure of low- noise amplifiers is of 1-2 dB, depending 
on technology, on frequency, and on bandwidth. The effect of bandwidth 
shows up clearly in microwave amplifiers, where the active devices have 
low noise temperature and capacitive input. The loss of the impedance- 
matching network necessary to match the resistive amplifier input (50 fi) 
to the capacitive transistor gate turns into increased noise figure. Needless 
to say, larger bandwidth design turns into higher input loss, thus into higher 
noise figure. 

The definitions of and F implicitly assume that the amplifier noise 
is a random process added to the useful signal, and not correlated to it. In 
the presence of a sinusoidal carrier of power Pq , the phase noise is 

FkT 

S^{f) = bo = — - — constant. (1.39) 

1.4.2 Flicker noise 

It has been experimentally observed |HWB68[ IWF.)97[ IHHWW();-i| that 
phase flickering of different amplifier types falls in a relatively narrow range, 
and that for a given amplifier the phase flickering is about independent of 
the carrier power. Consequently the flicker noise of m cascaded ampliflers of 
the same type is about m times the noise of one amplifier. This is radically 
different from the case white noise, where the noise of a stage referred to 
the input is divided by the gain of all the preceding stages. 

Table E21 shows the typical phase fiickering of commercial amplifiers. 
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Table 1.2: Typical phase flickering of amplifiers. 



RATE MOS bipolar HF 

microwave HF /UHF 



fair 

good 

best 



-100 
-110 
-120 



-120 
-130 
-140 



dBradVHz 



S(f) 



near-dc 
flicker 



no carrier 



S(f) 



parametric 
up-conversion 




Figure 1.12: Parametric up-conversion of near-dc flicker in amplifiers. 



The mechanism that originates phase flickering is a low-frequency (close 
to dc) random process with spectrum of the flicker type that modulates 
the carrier (Fig. ITTT^ . This mechanism is often called parametric noise^ 
because the near-dc flickering modulates a parameter of the device high- 
frequency model. 

A simplified nonlinear model provides insight in the counterintuitive 
behavior of phase noise. Let us consider a signal x{t) = Vq cos{ujot) + n(t) 
inside a device, where the sinusoid is the input signal, and n{t) is the near- 
dc flickering of the dc bias. Additionally, let us assume that the device 
is slightly nonlinear, and that the nonlinearity can be expanded as the 
2nd-degree polynomial P{x) = + aix + a2X^ . Feeding x{t) into the 
polynomial, we get a carrier term aiVo cos{uj{)t) , plus a close-in noise term 

■^Of course, the term parametric noise is more general than the phase flickering. 
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a2Von-{t) cos{uJt). Inspecting on the radio-frequency spectrum around loq, we 
find that the power of the noise sidebands is proportional to Vq, thus to the 
power of the input signal. This is in agreement with the simple observation 
that there can not be close-in noise in the absence of a carrier. As a 
consequence, the fractional amplitude fluctuation, i.e., the close-in noise 
divided by the carrier, is a{t) = ^n(i), independent of the carrier power. 
Similarly, a 2nd-degree nonlinear model that involves a variable reactance 
mechanism (varactor effect in transistors and MOS) leads to close-in phase 
noise independent of the carrier power. 

Phase flickering depends on the physical size of the amplifier active 
region. This can be proved through a gedankenexperiment in which we split 
the input signal into n equal branches, amplify and recombine. The power 
gain is that of a single branch. Conversely, the flicker noise is (6_i)tot = 
^(6_i)branch because the branch amplifiers are independent and the phase 
flickering of each is not affected by having reduced the power by n. If we join 
the n amplifiers in a single one that employs an active region (base in the 
case of bipolar transistor, or channel in the case of field-effect transistors) n 
times larger, the phase noise reduction is kept. The additional hypothesis 
is required, that the near-dc flickering takes place at microscopic scale, for 
there is no correlation between the different regions of the n-volume device. 
This hypothesis is consistent with the two most accredited models for the 
flicker noise |Hoo69[ ITvEWHT] . 

Finally, phase flickering is related to the amplifier gain. This is a side 
effect of the number of stages needed for a given gain, rather than a gain 
effect in a single stage. 

1.4.3 Other noise types 

Noise phenomena with a slope significantly steeper than —1 in the spec- 
trum, say /~^, do not exist in amplifiers. If such phenomenon was present, 
the delay of the amplifier would diverge rapidly. This never happens in 
practice. This general statement does not exclude some bumps in the spec- 
trum, for example due to the environment temperature, which in some 
frequency range yield a spectrum locally steeper than f~^. 
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Figure 1.13: Typical phase noise of an amplifier. 



1.4.4 Phase noise spectrum 

The total phase noise spectrum results from adding the white and the flicker 
noise spectra, as in Fig. 11.131 This relies on the assumption that white and 
flicker phenomena are independent, which is true for actual amplifiers. 

It is important to understand that bo (white) is proportional to the 
inverse of the carrier power Pq, while 6_i (flicker) is about independent of 
Pq. The corner frequency fc depends on the input power. The belief that 
fc is a noise parameter of the amplifier is a common mistake. 

1.4.5 Noise-corrected amplifiers 

The amplifier flicker noise is of paramount importance to the oscillator 
frequency stability. Reducing the flicker, even if this is done at expense of 
higher white noise, results in improved oscillator stability. A new generation 
of oscillators make use of a noise-corrected amplifier in the loop. This 
technique is based on a bridge scheme that takes the difference between 
the input and output of the amplifier, compensating for gain and delay. 
Balancing the bridge, the differential signal is the amplifier distortion and 
noise. The latter is amplified, detected, and used to compensate for the 
amplifier noise in closed loop. 

In a noise-corrected amplifier the input is split into two branches, that is. 
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Figure 1.14: With a flickering amplifier, the Leeson effect yields two types 
of spectrum, type 1 for fc < /l and type 2 for fc > fi- A noise- free 
resonator is assumed. 



the main amplifier and the noise-correction amplifier. A directional coupler 
may be necessary, which introduce a loss of 3 dB, plus the dissipative loss. 
Such loss turns into an increase in noise figure of the same amount. If the 
noise figure of the internal amplifier is of 1-2 dB and the total loss of the 
coupler is of 4 dB, the noise figure of the corrected amplifier is of 5-6 dB. 
On the other hand, the correction results in a reduction of 20-40 dB in the 
phase flickering. 

1.5 The phase noise spectrum of real oscillators 

For a given amplifier, the phase noise (Fig. II. ll?]) is white at high /, and of 
the flicker type below the cutoff frequency fc that results from the carrier 
power Pq. When such amplifier is inserted in an oscillator, it interacts with 
the resonator in the way shown in Fig. 11.141 Two basic types of interaction 
are possible. 

Type 1 is the most frequently encountered. It is typical of microwave os- 
cillator and high-frequency (> 100 MHz) piezoelectric oscillators, in which 
/l is made high by the high resonant frequency and by the low merit factor 
Q. Looking at Fig. ll.L^ from right to left, the amplifier phase noise is white 
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Figure 1.15: Effect of the resonator frequency fluctuations on the oscillator 
noise. 



and the Leeson effect originates white frequency noise (6_2/~^). At lower 
frequencies the amplifier phase noise turns into flicker, hence the oscillator 
noise turns into frequency flickering (6_3/~^). No flicker is present in the 
output phase spectrum. 

Type 2 is found in low-frequency (5-10 MHz) high-stability quartz os- 
cillator, in which the merit factor may exceed 10^. Looking at Fig. 11.131 
from right to left, the amplifier phase noise turns from white to flicker at 
f = fc- Accordingly, phase flickering is visible at the oscillator 

output. At lower frequencies the Leeson effect takes place, hence the os- 
cillator noise turns into frequency flickering (6_3/~^). No white frequency 
noise (6_2/~^) is present in the output phase spectrum. 

By inspection on Fig. ll.ltH there can be either the or the noise 
types, not both. 

1.5.1 The effect of the resonator noise 

The dissipative loss of the resonator, inherently, originates white noise. Yet, 
the noise phenomena most relevant to the oscillator stability are the flicker 
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and the random walk of the resonant frequency i^q. Thus, the spectrum 
Sy{f) of the fractional frequency fluctuation y = Av/iyQ shows a term 
for the frequency flicker, and /i_2/~^ for the frequency random 
walk. The relationship^ between S^{f) and Sy{f) is 

Syif) = 4 SM) ■ (1-40) 

^0 

Accordingly, the term h^if^^ of the resonator fluctuation yields a term 
proportional to in the phase noise, and the term /i_2/~^ yields a term 
/~^. The resonator fluctuation is independent of the amplifier noise, for it 
adds to the oscillator noise. 

Figure II. 151 shows the two basic spectra of Fig. 11.141 after introducing 
the resonator frequency fluctuation. The resonator fluctuation may hide the 
amplifier corner frequency fc, the Leeson effect, or both. Three behavior 
types deserve attention. 

1. The resonator noise hides fc but not /l- This is found in type-1 
spectra. The cross point of the resonator noise and the noise 
due to the Leeson effect has the same graphical signature of fc, but 
is not. 

2. The resonator noise hides fc and fi,. This is found in type-1 and 
type-2 spectra. Only one corner is visible on the plot, where the f~^ 
resonator noise crosses the f^ amplifier noise. This behavior is found 
for example in VHF quartz oscillators. 

3. The resonator noise hides fi, but not fc- This is expected^ in type-2 
spectra. 



1.5.2 The effect of the output buffer 

It is interesting to analyze the case of an oscillator with ideal resonator and 
noisy amplifier, followed by a noisy output buffer (Fig. Il.l(i|) . Of course 

^Start from the fractional frequency fluctuation y — j^^-^^- The time-domain 
derivative maps into a multiplication by u^, thus by (27r/) , in the spectrum. 
^This case needs some more theoretical analysis. 
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Figure 1.16: Effect of the output buffer on the oscillator noise. 



the output buffer is independent of the oscillator, for the two noise spectra 
barely add. 

The white noise comes from the additive white noise referred to the 
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carrier power Pq. The noise FkTo of the sustaining amphfier is amplified 
at the input of the buffer. Therefore, the white phase noise of the buffer 
ends up to be neghgible in virtually all practical cases. On the other hand, 
the flicker phase noise of amplifiers is about independent of the carrier 
power, hence the phase noise of the buffer is not expected to be negligible. 
Accounting for the buffer noise, the spectra of Fig. 11.141 are to be modified 
as shown in Fig. 11.161 

In type lA and 2A spectra, the sustaining amplifier and the buffer have 
similar flicker characteristics. In type lA (/^ > /c), the phase flickering of 
the buffer is hidden by the Leeson effect. The insertion of the output buffer 
lets the spectrum is substantially unchanged. 

In type 2A [fi < fc), the phase flickering of the buffer adds to the phase 
flickering of the oscillator. As a consequence, the corner point at which the 
noise turns into is pushed towards lower frequencies. This corner 
point can be easily mistaken for the Leeson frequency because it has the 
same graphical signature. 

In type IB spectrum, the sustaining amplifier exhibits low phase flick- 
ering, significantly lower than that of the buffer. As a result, fc is a low 
frequency. This is the case of some low- noise microwave oscillators, thus we 
assume fi > fc even in the case of large merit-factor oscillators. A simple 
amplifier is used as the buffer. If the merit factor of the resonator is large 
enough (low fi), the f~^ phase noise of the buffer shows up, hiding the f^ 
to slope transition characteristic of the Leeson effect. Yet, the Leeson 
frequency can be estimated extrapolating the segment still visible. 

If we assume that /l < fc, the case of the noise corrected amplifier 
and simple buffer is similar to the type 2A, but for a larger buffer noise. 
Consequently, the f~^ — > corner slides leftwards, farther from the true 
Leeson frequency. 



Chapter 2 



Oscillator hacking 



The combined knowledge of the oscillator noise theory, of general physics, 
and of electronic technology, enables to understand the inside of an oscil- 
lator from the data sheet, and to guess some relevant internal parameters 
like Pq, Q, /l, amplifier 1/f noise, etc. We describe the guidelines of this 
process and show some examples. The need of guessing the internal tech- 
nology is a source of difficulties and of inconsistencies which may make 
the interpretation only partially reliable. Coping with this is a part of the 
message addressed to the reader. 

Inspection on the data sheet. The first step consists of reading care- 
fully the data sheet focusing on the resonator and on the amplifier tech- 
nology, and bringing up to the mind as many related facts as possible. For 
example, a 5 MHz quartz can have a merit factor in excess of 10^, but it 
must be driven at very low power, say 10-20 /iW, for best long-term sta- 
bility. The merit factor of a dielectric resonator can be of 1000 or more, 
depending on size and frequency. So on for the other resonator types, and 
for the amplifiers. Similar oscillators encountered in the past may have a 
similar spectrum, or be surprisingly different. 

Parametric estimation of the spectrum. This part of the process 
consists of matching the phase noise spectrum with the polynomial S^{f) = 
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Yl^=-ibif^ in order to identify the coefficients bi. A term bif^ on a log-log 
plot appears as a straight line 



noise type 


term 




slope 


white phase 


^0 
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The actual spectra are of the form 5'<^(/) = Yl^=-4: ^iP + ^jif)^ where 
the terms Sj{f) account for the residuals of the mains (50 Hz or 60 Hz 
and multiples), for bumps due to feedback and for other stray phenomena. 
Figures 12.11 to 12.91 provide some examples of actual phase noise spectra. 
The mathematical process of matching the spectrum to a model is called 
parametric estimation |PW98[ IJW68j . Some a-priori knowledge of the na- 
ture of the stray signals may be necessary to match the complete model 
Yl^i=-i biP + ^j ^jif) to the observed spectrum. Although (almost) only in 
the power-law coefficients bi are relevant in the end, the Sj(/) are essential 
in that they reduce the bias and residuals of the estimation. 

Whereas computers provide accuracy, a general parametric estimator is 
not easy to implement. Conversely, the human eyes do well in filtering out 
the stray signals and getting a good straight-line (polynomial) approxima- 
tion. The inspection on a log-log plot by sliding old-fashion squares and 
rulers proves to be surprisingly useful. Spectra are often shown as £(/), 
don't forget that £(/) = ^Sip{f). Pick up the reference slope /* using the 
largest possible area on the coordinate frame, slide the square until its side 
coincides to the corresponding portion of the spectrum (6i/*), and get the 
coefficients bi. Proceed from right to left, thus from bo to 6-4. 

Generally, at the corner between two straight lines the true spectrum 
is 3 dB above the corner point. This is due to either one of the following 
reasons. In the first case, the difference in slope is 1 at the corner point. 
When this occurs, there are two independent random processes whose spec- 
trum takes the same value (6i/* = bi^if^^^) at the corner point. In the 
second case the difference in slope is of 2, for at the corner point it holds 



February 2, 2008 



E. Rubiola The Leeson effect 



29 



that bif^ = 5j+2/*^^. This occurs when a single noise process is filtered, 
due to the Leeson effect. The factor 2 (3 dB) at the corner point results 
from a single real zero of the complex transfer function. 

In some cases the difference between the spectrum and the straight-line 
approximation at a corner is not of 3 dB. When this occurs, one should 
work two estimations, one based on the straight-line fitting, and the other 
based on the 3 dB difference between straight lines and true spectrum at 
the corner frequencies. The best estimate is a weighted average of the two 
worked-out spectra. Physical judgment should be used to assign unequal 
weights. 

Interpretation. This part of the process starts from the identification of 
the spectrum type, among those analyzed in the previous chapter. Then we 
get into the learning process, in which each oscillator is a unique case. As a 
general rule, one should proceed from the right-hand side of the spectrum 
to left, thus from white phase noise to frequency flicker or to random walk. 

Starting from the white phase noise, we evaluate the power Pq at the 
input of the sustaining amplifier using S^{f) = bo = FhTo/Po [Eq. (|1.39() ]. 
Thus, 



One may admit a noise figure F = 1 dB for conventional amplifiers, and 
of F = 5 dB for noise-corrected amplifiers, due to the input power splitter. 
Thanks to the gain of the sustaining amplifier, the white noise of the output 
buffer can generally be neglected. 

The next step is to evaluate fc (flicker of the sustaining amplifier) and 
/l, in order of occurrence from right to left. It is then necessary to guess 
the oscillator sub-type (Fig. I1.16|) . A major difficulty is to understand 
whether the oscillator stability derives from the Leeson effect or from the 
resonator fluctuation. Inverting Eq. H1.28() . the Leeson frequency gives the 
merit factor 
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Figure 2.1: Phase noise of the 10 GHz DRO Miteq D210B. Courtesy of 
Miteq Inc. Interpretation, comments and mistakes are of the author. 



The corner frequency fc reveals the phase flickering of the amplifier 

(fc-l)ampli = bofc (2.3) 

If the the spectrum is of the type IB (noise-corrected sustaining amplifier), 
the 1/f noise is the phase noise of the output buffer. 

Finally, the Allan variance of the fractional frequency fluctuation 
(Tab. [TTT|) . i.e., the oscillator stability, is 



a,^(r) = ... + --#- + 21n(2)-# + — -^r + ... (2.4) 



2.1 Miteq DRO mod. D-210B 

Figure 123 shows the phase noise spectrum of the dielectric-resonator oscil- 
lator (DRO) Miteq D-210B, taken from the device data sheet. The plot is 
fitted by the polynomial Yl^=~3 ^if^, with 
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bo -146 dB 2.5x10"^^ radVHz 

6_i (not visible) 

b-2 -11 dB 7.9xl0~2 radVHz 

b-3 +37 dB 5.0x10^ rad^/Hz 

This indicates that the spectrum is of the type lA of Fig. 11.161 

One might be tempted to fit the spectrum with a smaller bo (say, —147 
dBrad^/Hz) and to add a term b-if~^ tangent to the curve at / ~ 2 MHz. 
In this case the spectrum would be of the type IB, which contains the signa- 
ture of the output buffer. We discard this alternate interpretation because 
the noise of the output buffer would be 6_i « 10"^ rad^/Hz (—80 dB), 
which is too high for a microwave amplifier (Table [T3|) . 
The spectrum gives the following indications. 

1. The coefficient bo derives from the amplifier noise FkT referred 
to the carrier power Pq at the input of the sustaining amplifier, 
that is, 6o = Assuming that the noise figure is F = 1 dB, 
thus FkT = 5.1x10-21 rad^/Hz (-173 dBradVHz), it follows that 
Po = ^ = 2 (-27 dBm). 

2. However arbitrary the assumption F = 1 dB may seem, it is repre- 
sentative of actual microwave amplifiers. Depending on bandwidth 
and technology, the noise figure of a "good" amplifier is between 0.5 
dB and 2 dB. In this range, we find Po between 1.8 /uW and 2.5 fiW. 

3. The spectrum changes slope from to /"^ at the Leeson frequency 
/l — 4.3 MHz. At this frequency, the asymptotic approximation is 
some 4 dB lower than the measured spectrum, instead of the expected 
3 dB. This discrepancy is tolerable. From fi ~ 4.3 MHz, it follows 
that Q = ~ 1160, quite plausible for a dielectric resonator. 

4. The white frequency coefficient is 6_2 = 7.9x10-^ rad^/Hz 
(-11 dBrad^/Hz). 



5. The corner point at which the slope changes from —2 to —3 is 70 kHz. 
This is the corner frequency fc of the amplifier, at which it holds that 
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(^-i)ampii/ ^ = (6o)ampii- Hence (6_i)ampii = 1.8x10 ^° rad^/Hz 
(-98 dBradVHz). 

6. The flicker frequency coefficient is 6_3 = 5x10^ rad^/Hz 
(+37 dBradVHz). 

7. Tlie wliite and flicker frequency noise, transformed into Allan variance 
(Table EH), is 

a2(r) = ^ + 21n(2) h.^ 

7.9x10-2 1 5x10^ 
+ 1.39X- 



2 X (1010)2 r (1010)2 ' 



thus 



^2(^)^ 4x10 '% g9^io-i7 
y r 

, , 2x10-11 q 

ay{T) ~ ^— + 8.3x 10"^ . 

V''" 

Finally, one should note that the oscillator flicker shows up in the 1-100 kHz 
region. Common sense suggests that temperature and other environmental 
fluctuations have no effect at this time scale, and that the flickering of 
the dielectric constant in the resonator will not exceed the amplifier noise. 
Consequently, in this region the oscillator flicker is due to the amplifier 
through the Leeson effect, rather than to the resonator. 



2.2 Poseidon DRO-10.4-FR (10.4 GHz) 

The Poseidon DRO-10.4-FR is another example of oscillator based on a 
dielectric resonator. Figure 12.21 shows the phase noise spectrum, from a 
preliminary data sheet. S^{f) is fitted by the polynomial Yl^=~3 ^iP^ with 
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Figure 2.2: Poseidon DRO 10.4-FR. The spectra are plotted using the 
preliminary data available on the manufacturer web site. Interpretation, 
comments and mistakes are of the author. 
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Once again, the spectrum is of the type lA fFig. ll.lUI) . typical of microwave 
oscillators. Yet, the discrepancy with respect to the theoretical model is 
larger than in the case of the Miteq oscillator. The spectrum gives the 
following indications. 

1. The spectrum results from the comparison of two DRO-10.4 oscilla- 
tors. In the absence of other indications, we believe that the numerical 
data refer to a single oscillator, after taking away 3 dB from the raw 
data. 

2. The white phase noise is 6o = 3.2xl0~i'^ radVHz (-165 dBradVHz). 
Thus Po = ^ - 160 (-8 dBm), assuming that F = 1 dB. 
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3. The Leeson frequency is Jl = 3.2 MHz. Accordingly, the merit factor 
is Q = 7JJ- — 1625, which is reasonable for a dielectric resonator. 

4. In a type-lA spectrum it holds that 6^2/"^ = bo at f = f^. Thus, the 
white frequency noise is 6_2 ^ 3.2x10"'^ radVHz (-35 dBradVHz). 

5. There is some discrepancy between the Leeson model and the true 
spectrum. In the region from 2 kHz to 200 kHz, the spectrum seems 
to be close to a line of slope /~^^^, rather than or /~^. At 
the present time this discrepancy, (up to 4 dB at / « 10 kHz) is 
unexplained. 

6. The corner frequency of the amplifier (i.e., the frequency at which 
the oscillator spectrum changes from to f~^) is fc = 9.3 kHz. 
Accordingly, the phase noise spectrum of the amplifier, on the 
left hand of / = fc, is (6-i)ampii = 60/c = 2.9x10-13 rad^/Hz 
(-125 dBradVHz). 

7. The amplifier flickering, 5 dB lower than the best in Table 11.21 is 
surprisingly low for a microwave amplifier contained a commercial 
product. Such a low noise could be obtained with SiGe technology, 
with a single-stage amplifier employing a large-volume transistor, or 
with some feedback or feedforward noise degeneration scheme. A 
noise degeneration scheme seems incompatible with the size of the 
packaged oscillator. Yet, nothing can be taken for sure on the basis 
of the available information. 

8. The flicker frequency coefficient is 6_3 = 2.5 rad^/Hz 
(+4 dBradVHz). 

9. The white and flicker frequency noise, transformed into Allan variance 
(Table [HI), is 
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2.5 1 5x10^ 

+ 1.39X- 



2 X (10.4x109)2 r ■ (10.4x109)2 



thus 



1.5x10-24 
r 

1.2x10 



cj^(r) ~ - — + 3.2x10"^° 



r 

12 



aJr) ~ + 1.8x10^^° . 

Vt 

Figure [221 also reports the phase noise spectrum of the DRO-10.4-XPL os- 
cihator, which is a different version of the same base design, intended for 
phase-locked loops. Below a cutoff frequency of about 70 kHz, this oscilla- 
tor is locked to an external reference, for the spectrum gives no additional 
information in this region. Nonetheless, it is to be noted that the spec- 
trum is proportional to f~^^'^, i.e., —25 dB/decade below the loop cutoff 
frequency of 10 kHz. This is the signature of a fractional-order control 
system, like that proposed in |(]M85j . 



2.3 Poseidon Shoebox (10 GHz sapphire res- 
onator) 

The Poseidon Shebox integrates a sapphire whispering gallery (WG) res- 
onator and an interferometric noise degeneration scheme. There results 
a low-noise oscillator intended for high short-term stability applications. 
Figure Eini shows the phase noise spectrum with a tentative interpretation. 
The spectrum seems to be of the type 2 fFig. II. ll)) . with Jl < fc- Qualita- 
tively, this is consistent with the fact that the WG resonator features high 
Q. Yet this interpretation suffers from three problems. 

1. From Q = we get Q f« 1.9x10^ at z/q = 10 GHz. This value is 
incompatible with the dielectric loss of the sapphire. For comparison, 
the typical merit factor of a 10 GHz WG resonator is in the range of 
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100 1000 fL-2.6kHz jq^qq 100000 



■-I Fourier frequency, Hz 

Figure 2.3: Poseidon Shoebox, first attempt to interpret the spectrum. 
Data are taken from tlie manufacturer data sheet. Interpretation, com- 
ments and mistakes are of the author. 



2x10^ 295 K room temperature 
3x10"^ 77 K hquidN 
5x10^ 4K hquidHe 

The dielectric loss of the sapphire is a reproducible function of tem- 
perature. Thus, the loss of the resonator depends on the space dis- 
tribution of the electric field, i.e., on the mode, in a narrow range. 
A merit factor of 1.9x10^ can not be obtained by moderate cooling 
(Peltier cells). Finally, the size and weigh (3 dm^ and 6.5 kg) in- 
dicate that the oscillator works at room temperature, at most with 
temperature stabilization. 

2. The phase flicker of the sustaining amplifier is with 6_i ~ 

10"^2 rad^/Hz. This is seen on the plot, -160 dBrad^/Hz at 
/ = 10 kHz with slope /~^ thus -120 dBrad^/Hz at 1 Hz. This 
value is too high for a sophisticated amplifier that makes use of the 
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Poseidon Shoebox 

10 GHz sapphire WG resonator 

noise correction 



T3 



J= - 




100 



1000 10000 
Fourier frequency, Hz 



Figure 2.4: Poseidon Shoebox, revised interpretation. Interpretation, com- 
ments and mistakes are of the author. 



interferometric noise correction technique. 

3. At fi, where b-3f~'^ = the true spectrum differs from the 

asymptotic approximation by 6 dB instead of 3 dB. This discrepancy 
is an additional reason to reject the interpretation of Fig. 12.31 

The above difficulties make us understand that the spectrum is of the type 
2A (Fig. 11.16(1 . in which the flicker noise of the output buffer shows up in 
the region around fi. There follows new interpretation, shown in Fig. 12.41 

1. As usual, we start from the white noise floor. On the Figure, we 
observe that bo = 1.3x10"^^ rad^Hz (-169 dBradVHz). This is 
ascribed to the sustaining amplifier. 

2. The sustaining amplifier makes use of an interferometric noise degen- 
eration circuit to reduce the fiicker noise. In this circuit there are two 
amplifiers, the first amplifies the input signal, and the second ampli- 
fies the noise of the first, at the output of a carrier suppression circuit. 
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The second amplifier corrects for the noise of the first by means of 
a feedback circuit. We guess a noise figure F = 5 dB, which results 
from the intrinsic loss of the power splitter at the input of the noise 
corrected amplifier (3 dB), from the resistive loss of the power split- 
ter and of the lines (1 dB), and from the noise figure of the second 
amplifier (1 dB). 

3. From bo = ^,we get Pq = 1 mW (0 dBm). 

4. The phase flickering of the output buffer shows up on the right- 
hand part of the spectrum, at 10^-10^ Hz. The noise coefficient is 
(6„i)buf ^ 10"^2 radVHz (-120 dBradVHz). The output buffer is 
a good microwave amplifier. 

5. After removing the buffer phase noise, the white frequency noise 
6_2/~^ of the oscillator is clearly identified. The decade centered 
at 4 kHz is used to find 6_2 = 7.9x10-^^ rad^/Hz (-81 dBrad^/Hz) 
on the plot. 

6. From fe_2/~^ = 6o) found on the plot, it follows that /l ~ 25 kHz. 
Hence, Q = — 2x10^. This is the typical value for a room- 
temperature WG sapphire resonator. 

7. It is seen on the plot that 6„3 ~ 7.9x 10"*^ rad^/Hz (-51 dBrad^/Hz). 

8. In the spectrum of the type 2A, the corner frequency fc of the sus- 
taining amplifier is the frequency at which the oscillator noise changes 
slope from to /~^. Thus, fc — 850 Hz. The flicker noise of the 
sustaining amplifier is (6-i)ampii/~^ with (6_i)ampii = &o/c- Thus, 
(^-i)ampii - 10-^^ radVHz (-140 dBrad^/Hz). 

9. The white and flicker frequency noise, transformed into Allan variance 
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(Table EH), is 



a^(r) = ^ + 21n(2) h^^ 

,2 
'0 



7.9x10"^ 1 7.9x10^6 
- + 1.39X- 



2 X (10x109)2 r (10x109)2 ' 



thus 



n, , 4x10-29 

y T 

ajr):^ ^-^^^"''' +3.3x10-^3 _ 
^/r 

10. The bump at / = 40 kHz is ascribed to the feedback control of the 
noise degeneration circuit. 

2.4 Oscilloquartz OCXO 8600 (5 MHz AT-cut 
BVA) 

Figure 12.51 shows the phase noise of the oven-controlled quartz oscillator 
Oscilloquartz OCXO 8600. This oscillator has been chosen as an example 
because of its outstanding stability in the 0.1-10 s region. The spectrum 
is of the type 2 {/l < fc), typical of low-frequency oscillators with high-Q 
resonators. The plot is fitted by the polynomial Yl^=-3^iP' with 

bo -155 dB 3.2x10-16 rad2/Hz 

b-i -131 dB 7.9x10-1^ rad2/Hz 

b-2 (not visible) 

6_3 -124 dB 4x10-13 rad2/Hz 



A preliminary interpretation the spectrum as follows. 
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f =250Hz 



Figure 2.5: Phase noise of the oven-controhed quartz oscillator Oscillo- 
quartz OCXO 8600, with a preliminary (wrong) interpretation. Courtesy 
of Oscilloquartz S.A. Interpretation, comments and mistakes are of the 
author. 



1. The white phase noise is 6o = 3.2x10-^^ rad^/Hz (-155 dBrad^/Hz). 
Thus 

FkT 

Po = — — ~ 16 (-18 dBm) , 

oq 

assuming that F = I dB. This power level is consistent with the 
general experience on 5~10 MHz oscillators, in which the power is 
kept low for best stability. 

2. /c ~ 250 Hz (from the plot), hence (6_i)ampii = 7.9 xlQ-^^. 

3. On the plot we observe that ~ 2.2 Hz (/^ in Fig. I2.5jl . and 6„3 ~ 
4x10-^3. 

4. g = hence Q ~ 1.1 xlO^. 
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(guessed) 



f =80Hz 



Figure 2.6: Phase noise of the oven-controlled quartz oscillator Oscillo- 
quartz OCXO 8600, with revised interpretation. Courtesy of Oscilloquartz 
S.A. Interpretation, comments and mistakes are of the author. 



The item 1 is correct because the effect of the output buffer on the white 
noise is divided by the gain of the sustaining amplifier, thus it is negligible. 
The white noise originates in the series resistance of the resonator and in 
the sustaining amplifier input. Conversely, the flicker noise of the output 
buffer is not negligible because it results from parametric modulation. Thus 
we change the interpretation (Fig. 12. 6 j) as follows. 

1. Section fl.5.21 indicates that the effect of the output buffer can not be 
neglected if phase flickering shows up in the phase-noise spectrum. 
We guess that there are two stages between the sustaining amplifier 
and the output and that the 1// noise of each stage is at least equal 
to that of the sustaining amplifier. Thus the 1/f noise of the sus- 
taining amplifier is not higher than 1/3 of the total noise. Let us 
take away 5 dB. In reality, the difference could be more then 5 dB 
because the output amplifier can be more complex, and because su- 
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perior technology can be used in the sustaining amplifier, rather than 
in the buffers. 

2. The new estimate of the flickering of the sustaining amplifier is 
(^-i)ampii = 2.5x10-1^ radVHz (-136 dBradVHz). This affects 
the corner frequency and in turn the Leeson frequency. The new 
values are 

fc ^80Rz 

and fi ~ 3.8 Hz (/^ in Fig. 12. 6|) . The latter is still provisional. 

3. Using fi = 3.8 Hz in /l = ^, we get Q ~ 6.6x10^ 

4. Q ~ 6.6x10^ is a too low value for a top-technology oscillator. The 
literature suggests that the product uqQ is between 10^^ and 2x10^^. 

5. Q is higher than value of 6.6x10^. This means that the Leeson effect 
is hidden, and that we are not able to calculate Q. Henceforth, our 
knowledge about Q relies only upon experience and on the literature. 

6. If we guess Q = 2.5x10^, reduced to Q = 1.8x10^ (loaded) by the 
dissipative effect of the amplifier input, we get 

= 1.4Hz 

This indicates that the oscillator frequency flicker 6_3 ~ 3.2x10^^ 
is due to frequency flickering in the resonator. The Leeson effect 
is hidden by the resonator instability and by the phase noise of the 
output buffer. Thus, the corner frequency f'l ~ 3.8 Hz of Fig. 12.61 is 
not the Leeson frequency. 

7. The flicker noise of the sustaining amplifier combined with the Leeson 
frequency gives the stability limit of the electronics of the oscillator. 
This is the solid line in Fig. 12. 6| some 10 dB below the phase 
noise. 
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i. The frequency flicker turned into Allan variance (Table [TTT|) . is 
cj^(r) = 21n(2) = 21n(2) = 1.39x- 



thus 



(5x106)2 

cj^(t) ~ 2.2x10"^^ o-2^(r) ~ 1.5x10"^^ . 



This is lower than the value cry(r) < 3x10"^^ for 0.2s < r < 30s. 
given in the specifications. On the other hand, specifications are 
conservative, while a sample can be better. 



2.5 FEMTO-ST prototype (10 MHz LD-cut 
quartz) 

The interest of the LD cut is the low isochronism defect — ^td, i.e., the 
low dependence of frequency on the drive power ._GB B02j . For comparison, 
the isochronism defect can be of 1.2 x 10^^^ / /xW (typical) for the popular 
SC-cut crystals, and as low as 10~^^ / //W for the LD cut. This oscillator is 
a laboratory prototype. It has been included in this Section otherwise de- 
voted to commercial oscillators because we know some internal parameters, 
for the analysis process is somewhat different. 

Figure 12.71 shows the phase noise spectrum. The spectrum is the total 
noise obtained by comparing two equal devices, and divided by two. The 
plot is fitted by the polynomial Yl^=-3 biP, with 

6o -147 dB 2xl0~^5 

6_i -130 dB 1x10-13 

b-2 (not visible) 

b-3 -116.6 dB 2.2x10-12 

We interpret the spectrum as follows. 

1. The white phase noise is bo = 2x10""^ radVHz (-147 dBradVHz). 
As we know the power at the amplifier input, Pq = 340 /uW, we 



radVHz 
rad^/Hz 

radVHz 
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Fourier frequency, Hz 



calculated 



Figure 2.7: Phase noise of a FEMTO-ST oscillator prototype based on a 
LD-cut resonator. Courtesy of the FEMTO-ST laboratory. Interpretation, 
comments and mistakes are mine. 



can calculate the noise figure F = ^ ~ 167 (+22 dB). This value 
is far too high for the noise figure of a HF-VHF amplifier, thus the 
observed floor can not be the noise floor bQ of the oscillator. Therefore 
we ascribe it to the white noise floor of the instrument, used at low 
input power. 

2. As 6o is hidden, the corner frequency fc of the amplifier is also hidden. 

3. The phase flickering is clearly visible on a span of one decade. At 
the corner frequencies, 5 Hz and 50 Hz, the difference between exper- 
imental data and asymptotic approximation is of 3 dB as expected. 
The residual noise of the instrument is not shown on the plot. Yet ex- 
perience suggests that the value (6_i)osc = 10~^^ is sufficiently larger 
than the noise of a HF-VHF mixer. This elements indicate that the 
value (6_i)osc = 10~^^ can be trusted. 
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4. From the plot, the frequency flicker is 5_3 = 2.2xl0~ . We use it 
to estimate the Leeson frequency, at which b-sf"^ = There 
results the value = 8.9 Hz. 

5. Let us assume that the sustaining amplifier contributes 1/3 of the 
total phase flickering, guessing that there are two buffer stages (Sec- 
tion ^212) based on the same technology of the sustaining amplifier. 
There follows that (6-i)ampii = 3.16x10"^^. 

6. The merit factor of the resonator in the circuit load conditions is 
known, Q = 1.15x10^. The calculated Leeson frequency is 

fL = ^ = 4.35 Hz . 

This value suffers from the uncertainty with which we guess the phase 
flickering of the sustaining amplifier as a fraction of the total phase 
flickering. Nonetheless, the buffer is necessary, which ensures that 
only a fraction of the phase flickering is due to the sustaining amplifier. 

7. The calculated Leeson frequency (4.35 Hz) is sufficiently lower than 
the apparent value (8.9 Hz) to provide evidence that the Leeson effect 
is not visible on the spectrum, and that the 1//^ line is the frequency 
flickering of the resonator. 

8. The frequency flicker turned into Allan variance f Table [TTT|) . is 

a2(T) =21n(2) h_i 

hence 

a^(r) ~ 3xl0~26 ayir) ~ l.JixlQ-'^^ . 

9. The total 1//^ noise results from the amplifier noise through the 
Leeson effect and from resonator frequency flickering 



, , , 6_3 2.2x10-^2 

= 21n 2 = 2x- — 

^ ' vl 10x106 2 



(^-3)osc = (^-3)lc + (ft-3)rcson 



(2.5) 
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Thus we can infer and the resonator stabihty 



a2(r) = 21n(2)( 




osc 




) 



The numerical result 



cj^(r) ~ 2.2x10 



,-26 



ay{T) ~ 1.5x10 



-13 



should be taken as an indication only, rather than as a measure, 
because it suffers from uncertainty enhancement of the difference, 
and from the lack of reliable measurement of the sustaining amplifier. 

Remark 2.1 This example provides a method for the independent mea- 
surement of the resonator stability. The accurate measurement of the am- 
plifier fiicker (6_i)ampii and of the merit factor Q in the appropriate load 
conditions turns into a reliable measure of /l and of (6_3)Le to be fed in 



2.6 Wenzel 501-04623 (100 MHz SC-cut quartz) 



The Wenzel 501-04623 oscillator is chosen as an example because of its 
outstanding low phase noise fioor, and because of its high carrier frequency 
(100 MHz). These two features are necessary for low noise after frequency 
multiplication. A table of the values of phase noise is given by the manu- 
facturer, instead of the complete spectrum. These values are used to plot 
the spectrum (Fig. 12. 8j) 



Eq.EB 



▲ 



bo 
b-i 
b-2 
b^3 



-173 dB 



5x10-1*^ 
(not visible) 



rad^/Hz 



-97 dB 



(not visible) 
2x10"^° 



rad^/Hz 



and to interpret it as follows. 
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Figure 2.8: Phase noise of the 100 MHz oscillator Wenzel 501-04623. Data 
are from the manufacturer web site. Interpretation, comments and mistakes 
are mine. 



1. The white phase noise is bo = 5x10"^'^ rad^/Hz (-173 dBrad^/Hz). 
Thus 

Pq = Eh^L ~ 1 niW (0 dBm) , 
&o 

assuming that F = \ dB. This relatively large power is necessary for 
low short-term noise, although it may be detrimental to the medium- 
term and long-term stability. 

2. The terms 6-2/"^ and b-if~^ are virtually absent, as there is no 
room for them in the region between the two specified points, 1 kHz 
and 10 kHz. A confirmation comes from the measurement of some 
samples. Consequently, the point at which the slope changes by 2 
(either — > or f~^) is not visible. This indicates that 
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the Leeson effect is hidden by the frequency flicker of the resonator. 

3. Of course, the corner point between the flicker of frequency and the 
white noise, 6_3/~^ = bo, which occurs at / = 3 kHz (/^ in Fig. I2.8j) . 
is not the Leeson frequency. 

4. As the Leeson effect is hidden, we are not able to calculate Q. Relying 
upon the literature, we can guess that a high-stability 100 MHz quartz 
resonator can have a merit factor of 1.2x10^, reduced to Q = SxlO"^ 
(loaded) by the power dissipation at the amplifier input. Accordingly, 
the Leeson frequency is 

F. = ^^625 Hz. 

5. The flicker noise of a state-of-the-art HF/VHF amplifier is 
(f'-i)ampii ~ 3.2x10^1^ radVHz (-135 dBrad^/Hz). This value, 
combined with the Leeson frequency, gives the stability limit of the 
electronics of the oscillator. This is the dashed line in Fig. 12. 8| 
some 10 dB below the phase noise. 

6. The frequency flicker turned into Allan variance (Table [TTT|) . is 

h 9x1 0""'^^ 

alir) = 21n(2) h_i = 21n(2) ^ = 1.39x 



thus 



2 ■ (100X106)2 ' 

(T^(t) ~ 2.8x10"^^ ayir) ~ 1.7x10"^^ . 



2.7 Oewaves Tidalwave (10 GHz OEO) 

The unique feature of the opto-electronic oscillator (OEO) is that the fre- 
quency reference is an optical delay line instead of a traditional resonator. 
The basic structure consists of a loop in which the microwave sinusoid 
modulates a laser beam. After traveling through a long optical fiber, the 
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Figure 2.9: Phase noise spectrum of the Oewaves Tidalwave photonic os- 
cihator Courtesy of Oewaves Inc. Interpretation, comments and mistakes 
are of the author. 



microwave signal is reconstructed by a high-speed photodetector, and fed 
back into the amplifier. Yet the actual structure may be more complex. 
Beside laboratory prototypes, Oewaves is up to now the only one manufac- 
turer of this type of oscillator. Studying the noise in oscillators, this is a 
challenging example because the general experience earned with traditional 
oscillators is of scarce usefulness. The phase noise spectrum (Fig; 12.9)1 is 
fitted by the polynomial Yl^=-3^if^y with 

bo -138 dB 1.6x10^^'^ radVHz 

b-i (not visible) 

b-2 (not visible) 

b-3 -21 dB 7.9xl0~3 radVHz 

b-4 +11 dB 1.26x10^ radVHz 

1. In a delay- line oscillator with a delay r we expect (Chapter a 
"clean" comb of spectral lines at /j = i/r, integer i > 1. This struc- 
ture is not present in the spectrum of Fig; 12.91 Instead, there is a 
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series of smaller bumps and spectral lines. This indicates that the 
Tidalwave is not a simple delay-line oscillator, and that it contains 
some additional circuits that reduce or almost eliminate the spectral 
lines at = /x/r, integer The remedy suggested by the literature 
is the dual-loop oscillator, in which a second delay line, significantly, 
is used to remove the largest peaks. One conference article |EM03j 
provides some interesting details, yet insufficient to identify clearly 
the inside of this specific oscillator. 

2. The right-hand side of the spectrum, in the white-noise region, an 
horizontal bo is only a poor approximation of the spectrum. This is 
inherent in the delay-line oscillator. Discarding the peaks, we take 
bo w 1.6x10-14 radVHz (-138 dBradVHz). 

3. The equivalent noise figure used to estimate the microwave power 
must be referred to the lowest-power point of the circuit, which is 
the output of the photodetector. Noise includes at least the intensity 
fluctuation of the laser, the shot noise, and the noise of the microwave 
amplifier. Let us guess that F = 6 dB. Consequently, the microwave 
power Po = at the phodetector output is Pq = 1 /^^W (—30 dBm). 

4. One may think that the transition f~"^ around / = 7 kHz is the 
signature of the frequency flicker of the reference, as in the Wenzel 
quartz oscillator. Yet an alternate interpretation is possible in this 
case because fc and /l might be close to one another and fall in this 
region. 

5. A corner frequency fc^7 kHz with bo ~ 1.6xW~^'^ rad^/Hz requires 
that 6_i = 10^1° rad^/Hz (-100 dBrad^/Hz). The latter value can 
be ascribed to the microwave amplifier (Tab. ITT^ . or to the optical 
system. 

6. The delay r is related to the Leeson frequency by fi = Thus 
fi^"^ kHz requires that r ~ 23 ^s, hence a fiber length / = ~ 4.7 
km, which is likely. 
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7. The — > transition at / ~ 7 kHz can be ascribed either to the 
fact that fc ~ /l, or to the flickering of the delay. The analysis of 
the spectrum is not sufficient to say more. 

8. The frequency flicker and random walk, turned into Allan variance 
(Table inj), is 

C72(r) = 21n(2) h_i + — h-2T 



thus 



(1010)2 6 (IQiO) 



cr^(r) ~ 1.1x10-22 + 8.3x10"^^ T 
o-y(r) ~ 1.05x10-11 + 9x 10-^ ^/T 
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Chapter 3 



Phase noise and linear feedback theory 



Figure proposes a model for the oscillator loop in the domain of the 
Laplace transforms. The amplifier of gain A (constant) and feedback path 
/5(s) are the blocks we are familiar with. The signal Vi{s) at the input of the 
summer S allows initial conditions and noise to be introduced in the loop. 
Interestingly, if Vi{s) is a driving signal, Fig. IH. II models for the injection- 
locked oscillator. Figure 13.11 differs form the classical control theory in 
that the at the feedback input of S positive, and in that the output is 
taken at the amplifier input instead of at the output. This choice simplifies 
some equations without affecting the phase noise. The elementary feedback 
theory tells us that the transfer function is 

H{s) = ^ detoi H{s) (3.1) 

is 

From this standpoint, an oscillator is a system that has (at least) a pair 
of imaginary conjugate poles, excited with suitable initial conditions. One 
may also interpret frequency stability as the stability of the poles. We 
first study the properties of Eq. (|3.2j) when /3(s) is the transfer function 
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Figure 3.1: Oscillator loop. 



of a simple resonator, ruled by a second-order differential equation. Then, 
we extend the analysis to the phase space in the Laplace domain, which 
describes phase noise. 



3.1 Oscillator and Laplace transforms 

The normalized transfer function of a resonator (Appendix^ is 
«r> -^0 resonator h = 

''''"=9 poles h:yjy 

By replacing /3(s) in Eq. (|3.2jl . we get 

The above H{s) = ^||y is a second-order rational function with real coeffi- 
cients. It has two poles, either real or complex conjugates, depending on the 
gain A. The root locus, shown in Figure has the following properties. 

1. The poles are imaginary conjugates for A = 1. 

2. The poles are complex conjugates for 1 — ^ < A < 1 + ^, with 
(Tp = — real elsewhere. 
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A - Phase transfer function H(s) 



B - Root locus of the denominator of H(s) 



J CO" 

constant zeros poles function of A 



increasing A 



60 

u 
> 



JCO 



A=l 
oscillation 



A=l 




A=l 



Figure 3.2: A: Noise transfer function H(s) [Eq. (|3.4j) ]. and B: root locus 
of the denominator of H{s), as a function of the gain A. 



3. The complex conjugate poles are on a circle of radius ujq centered in 
the origin. 

The proof is given underneath. 

Property The poles are imaginary conjugates for A = 1. 

Proof. We first rewrite the denominator T>{s) in a more convenient form 

V{s) = {s-s,){s-s;)-A'^s 

= {s-ap- jujp){s -crp+ ju)p) - A^s 

= {S- (Tp)+UJp - A — S 

= — 2apS + dp + + 2AapS 
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and finally 

V{s) = s"^ + 2ap{A-l)s + io^ (3.5) 

To obtain the above, we have used the properties + Up = ujq, and ap = 
— The poles of H{s) are the solutions of 'D{s) = 0, thus 

si, S2 = -ap{A - 1) ± ^a2(^-l)2-a;2 (3.6) 



Condition for stationary oscillation. Stationary oscillation requires that 
the poles are imaginary conjugates 

A = 1 which yields s = ijcjo 

This is obtained by replacing A = 1 in Eq. 13.51 The discriminator 

AV = a'p{A-lf-u;l 

reduces to — Wg, and the real part —ap{A—l) of the solutions si, S2 vanishes. 

▲ 

Property The poles are complex conjugates for 1 — ^<A<1 + ^, with 
ap = real elsewhere. 

Proof. The poles of ()3.4|) are complex conjugates when the discriminator 
AV < 0, thus 

a^p{A-lf-u;i<0. 
We first solve AP = 0, which yields 



Ai, A2 = l± 



ap 



As the coefficient of A in AT> is positive, we conclude that 
AV <0 for ^1 < A < ^2 . 
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The poles of H(s) are real coincident when AD = 0. We replace A = Ai 
and A = A2 inV 

+ 2crp ^1 ± ^ - 1^ + u;2 ^ . 

Hence the poles are real coincident for + 2ujqs + = 0, thus 

S3, S4 = ibwo • 

▲ 

Property The complex conjugate poles are on a circle of radius ujq centered 
in the origin. 

Proof. The complex conjugate solutions of P = are 



s = -a,{A - 1) ± jy u;2 - ajiA - 1)2 . (3.7) 

This is an alternate form of Eq. (|3.6j) . which makes the square root real 
when AD < 0. The square distance of the poles from the origin is 

= [R{s)f + [-s{s)f 

= al{A-lf+u:l-al{A-lf , 

which simplifies as 

i?2 = ojI- (3.8) 

A 

Remark 3.1 That the poles of H{s) are on a circle centered in the origin 
has an important consequence in metrology. In the real world the gain 
fluctuates around the value A = \. For small fluctuations of A, the poles 
fluctuate perpendicularly to the imaginary axis. The effect on the oscilla- 
tion frequency is of the second order only. ▲ 
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phase impulse in 
the voltage domain 



cos[Q)t+5 (t)] 



P(t) 



resonator 



phase impulse in 
the phase domain 



5(t) 



b(t) 



resonator 



phase impulse in cos[£Ot+U(t)] 
the voltage domain 



resonator 



Pt,(t) 



I 
I 

o 



phase step in U(t) 
the phase domain 



resonator 



b,,(t)= b(t)dt 



Figure 3.3: Resonator response to a phase impulse S{t). 



Remark 3.2 The exact condition A = 1 can not be ensured without a 
gain control mechanism. The latter can be interpreted as a control that 
stabilizes the oscillator poles onto the imaginary axis. 

The gain control can be a feedback that sets A for the output voltage 
to be constant. This approach was followed in early times of electronics 
in the Wien bridge oscillator |Pac39j . Yet amplifier saturation proved to 
be an effective amplitude control, even in ultra-stable oscillators. When 
the amplifier saturates, the output signal is clipped more or less smoothly. 
Power leakage from the fundamental to the harmonics reduces the gain, 
and in turn stabilizes the output amplitude. The resonator prevents the 
harmonics from being fed back to the amplifier input. ▲ 

3.2 Resonator in the phase space 

We analyze the resonator phase response b{t), and its Laplace transform 
B(s), in quasi-stationary conditions. The resonator is driven with a sinu- 
soidal signal at the frequency lo, which can be the natural frequency loq or 
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any other frequency. In the time domain, the phase transfer function b{t) 
is the phase of the resonator response to a Dirac 6{t) function in the phase 
of the input. More precisely, b{t) is defined as fohows (Fig. \'A.'A\i 

Vi{t) = cos[a;ot + 6{t)] input with phase impulse 

Vo{t) = cos[u;o* + output response, defines b{t) . 

We find b{t) from the response bu{t) to the Heaviside function (unit step) 



U(t) = / 6(t)dt = I Heaviside (3.9) 

J-oo 1 t>0 



using the property of linear systems that the impulse response is the deriva- 
tive of the step response bif{t) 

b{t) = ^bu{t) . (3.10) 

We calculate the response in small signal conditions, using a phase step 
KU{t), K ^ 0, so that some functions can be linearized. This is physically 
correct because phase noise in actual oscillators is a small signal. 

The method (Fig. 13.4(1 consists of breaking the input sinusoid Vi{t) at 
the time t = 

Vi{t) = Vi{t) + u-'(t) Vo{t) = v'^{t) + vl{t), output, for t > 
v[{t) = Vi{t)U{—t) ^oi't) response to Wj'(t), for t > 

v'^{t) = Vi{t)U{t) v"{t) response to v'-{t), for t > , 

so that the phase step KU{t) can be introduced in the phase of v'-{t). The 
terms v'g{t) and v"{t) are, respectively, the switch-on and the switch-off 
transient response for t > 0. 

3.2.1 Input signal tuned at the exact natural frequency 

Property The impulse response of the resonator in the phase space is 

bit) = -e-*/- (3.11) 
^ ' r response ^ ' 
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17 
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Figure 3.4: A sinusoid with a phase step can be decomposed into a sinusoid 
switched off, plus a shifted sinusoid switched on. The step response is the 
linear superposition of the two responses. 



and its Laplace transform is 
Bis) = 



1 



ST+1 



transfer 
function 



(3.12) 



Proof. Let 

Vi{t) = cos [ujot + kU (t)] 

= cos{ujQt)U {—t) + cos{ujot + K)U{t) 

the input signal. The resonator response v'^{t) to the sinusoid switched off 
is the exponentially decaying sinusoid 

v'^{t) = cos(u;ot)e"*/^ t > , 

where r = ^ is the resonator relaxation time (Appendix^. Similarly, 
the response v'^{t) to the switched-on sinusoid is the exponentially growing 
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sinusoid 

v'^{t) = cos{ujot + k)[1 - e-^/^] t>0 

The total output signal is 

Vo{t) =v'^it) + v'^{t) t>0 

= cos{LOot)e~^ + cos(woi) cos(k) — sin(u;o*) sin(K 



1 - e~ 



= cos(a;oi) e T + cos(k) — cos(K)e ^ — sin(a;oi) sin(K) 1 — e 
For K ^ we replace cos(k) = 1 and sin(K) = k, thus 

Vo{t) = cos{ujQt) — K.sm{LjQt)[l — e"/"^] 
Deleting ujQt, the above can be seen as a slowly varying phasor 

= l+jK[l-e-*/^] K<.1. 



<S^V (t)\ 

The angle arctan ^^{^(t)} ' iioi^™alized on k, is the step response 

step response 



hu{t) = 1 - e-*/^ 
Taking the derivative, we obtain the impulse response 



b{t) 



(3.13) 



which is Eq. ()3.11() . The Laplace transform [Eq. ()3.12() ] is found in common 
mathematical tables. The proof is omitted. ▲ 

3.2.2 Detuned input signal 

Property The step response of the resonator in the phase space is 
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o 



step response b_u(t) 


tau=1 


Omega=2 / 






Omega = 0, 0.5, 1, 1.5, 2 


// / // Omega=0 






time 



Figure 3.5: Step response bu(t) evaluated for r = 1. 



and the impulse response impulse is 



b{t) = n sin(Ot) + - cos(Ot) 
r 



-t/r impu lse r esponse 
(Fig. m^) 



where 



= cji — Wo def. of . 



(3.15) 



(3.16) 



is the frequency offset (or detuning), i.e., the difference between the oper- 
ating frequency uJi and the natural frequency ljOq of the resonator. 

Proof. The resonator asymptotic response {t oo) to the detuned input 
sinusoid Vi{t) = cos(cJit), uji / ojq, is Vo{t) = cos(u;it + (/>), with Voo and 
4) constants (cfr. Appendix ^ d. I1U5|) . Without loss of generality, we can 
scale and shift the input signal 



Vi{t) 



1 



cos(u;it — (f) , 
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Figure 3.6: Impulse response b{t) evaluated for r = 1. 



in order to get the asymptotic response of the form 



Vo{t) = cos(a;it) t — > oo . 

We use the method of the small phase step KU{t). Yet we need some 

preliminary results before introducing the step. 

Input signal switched off at the time t = 0. The input signal 

<{t) = ^cos{LOit-<p)U{-t) , 

switched off at t = by U{—t), yields the output 

v'^{t) = cos(woi)e"*/^ t>0 . (3.17) 
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Input switched on at the time t = 0. It is known from the differential 
equation theory that, in the presence of the input signal 

v'/{t) = -Lcos{u;it-4>)U{t) 

switched on at f = by U{t), the output is of the form 

v'^{t) = acos{iOot)e~T + bsm{LOQt)e~T +ccos(u;it) + dsm{ujit) t> 

with a, b, c, and d constants to be determined. As for t ^ oo it holds 
that Vo(t) = cos(t(Jit), we conclude that c = 1 and d = 0. Then, a and 
b are found by assessing the continuity of the output at t = when the 
input signal is Vi{t) = y— cos{uJit — (j)). The continuity condition is met for 
a = —1 and 5 = 0. 

In summary, it holds the following input-output relationship 

. / , . <W =^cos(a;it-0)t/(t) 

m/out pair ^ , ^ (3.18) 

v'l,{t) = [- cos(a;oi)e"*/^ + cos(a;it)] U{t) 

and, as an obvious extension, 

. , , . =^^H^it-ct>)u{t) 

m/out pair i"^ , ^ (3.19) 

v'^{t) = [-sin(a;ot)e-*/" + sin(L^it)] C/(t) 

Phase step at t = 0. Let us define the input signal Vi{t) = v'^it) + v"(t) 
v[{t) = ^cosicvit- <P)Ui-t) 

v'l{t) = ^ cos(a;it - 4> + K)U{t) k < 1 . 

The input v'l{t) can be rewritten as 

v'l{t) = p3- [cos(a;it — (j)) cos k — sin(a;it — (j)) sin k] U{t) 

= ^ [cos(Lo'it — (j)) — K sin(ajit — (j))\U{t) k ^1 . 
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Using Eq. (|:-{.17j) and the pairs (|H.18j) and (jS^J, we get the output signal 



= cos(u;oi)e~*^^ 



+ [— cos{ojot)e ^'"^ + COs(Lo'lt)] 
+ K[sin(wot)e~*/'^ — sin(a;it)] 



t > 

response to v'^ 
response to v" , 1st part 
response to v" , 2nd part 



hence 



Vo{t) = cos{LOit) — Ksm{u>it) + Kcos{u>ot)e ^^'^ t > . 



(3.20) 



Let us define the detuning frequency O = cji — Accordingly, it holds that 
uji^ = — O, hence sin(a;oO = sin(wit — fit), and consequently 

sin(i^oi) = sin(6<Jit) cos(i7t) — cos(Lo'it) sin(r2t) . 

The output signal (|3.2m) . written in terms of and 0, is 

Vo{t) = cos(a;it) — Ksin(a;it)+ 

+ Ksin(a;it) cos(Ot)e~T — ft;cos(t(Jit) sin(Ot)e~T 

which simplifies as 



Vo(t) = cos(a;ii) 1 — Ksin(Ot)e 



Ksin(a;it) 1 — cos(ilt)e r . (3.21) 



Output phase response. Freezing the oscillation ujit, the output signal (|3.2H) 
turns into a slow-varying phasor 



Vo{t) = 1 + jK;[l - cos(m)e 



-t/rl 



K < 1. 



The angle arctan ^j^^^^, normalized on k, is the phase-step response 



buit) = 1 — cos(r2t)e 



step response 
(Fig.ESI) 



66 



E. Rubiola The Leeson effect 



February 2, 2008 



Using the property b{t) = ^bu{t) [Eq. we find the impulse response 



b{t) = n sm{nt) + - cos(Ot) 
L r 



-t/r imp ulse response 

(Fig.Eni) 



Remark 3.3 In some unusual cases, a resonator may need a pump signal, 
for the measurement of the resonator parameter is difficult or impossible 
if the resonator is not in the oscillator loop. This occurs in the domain 
of photonics with the coupled opto-electronic oscillator (COEO) |YM97j . 
In such case the phase step introduced in the oscillator loop proved to be 
a useful method for the measurement of the resonator parameters. A mi- 
crowave phase modulator (varactor) was inserted in series to the resonator 
and driven with a square-wave generator at a low frequency. A 



Property The phase transfer function of the resonator in the presence of 
a detuned driving signal is 

, , _ 1 s + ^ + Q'^t transfer function , , 

^^')-r (s + ^-jn) (s + ^+jn) (Fzg.\rj\and\m ' ^ ' 



Proof. The Laplace transform B{s) = is found using the Euler 

identities 



cos{nt) = i (e^™ + e-J™) 
sm{nt) = ^ (e^™ - e-^™) 

and the properties of the Laplace transform 



s + l/r 

£{e'''f{t)}=F{s-a) 
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resonator 
tuned at m. 



A - resonator duned 

B(s) 



B - resonator detuned 



J CO 



resonator B(s) 
detuned by n«a)^, A 



J CO 



X 



-™l(1+£12t2) 

Figure 3.7: Phase noise transfer function B(s) on the Laplace plane. 



Thus 

B{s) = L 
and finally 

B{s) = - 



1 

i2 



r 2 



transfer function 



r (s + i - (s + 1 + jVi) (Fig.O and El) 



Remark 3.4 The phase noise bandwidth of the resonator increases when 
the resonator is detuned (Fig. 13. 8 j) . This is related to the following two 
interesting facts. 



arg[/3(ja;)] 

Hence, off the exact resonance, the 



1. At 6<j / the slope 



is lower than at (Jq (Fig. IA.1|) . 
resonator becomes less dispersive. 



2. The step response is faster when the resonator is detuned, as it is seen 
in Fig. 13.51 



▲ 
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Figure 3.8: Phase noise transfer function \B{juj)\'^ [Eq. (|3.22|) ] evaluated 
for T = 1, thus = 1. We observe that the transformation t = instead 
of r = 1 (hence = 1), ^ = z^i — z/q, and to ^ f, lets the plot 
unchanged. In other words, if the resonator Leeson frequency is the = 1 
and one reads / instead of w on the horizontal axis, the curves are the phase 
response \B{jf)\'^ of the resonator tuned, detuned by ^Jl, by Jl, etc. 



3.3 Another derivation of the resonator phase re- 
sponse 

The low-frequency phase transfer function of the resonator is 
1 



2|^(M)| 



H{j{uj-uJi))+H{j{io + iOi)) 



(3.23) 



The derivation^ of the resonator phase response given here is more general 
than that shown in Section f3. 21 which holds only for resonator. 

^Thanks to Charles (Chuck) Greenhall, JPL, Pasadena, CA. 
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Proof. Let h{t) the voltage impulse response of the resonator and H(juj) 
its Fourier transform. Let Vi{t) = cos[iOit + 6 + ip{t)] the input voltage. The 
output voltage is 

Vo(t) = {h* Vi){t) convolution 
h[u)vi{u) du 

L J — oo 

Linearize for (/? <C 1 

Vo{t) = 3f?|e^'('"^*+^) r h{u)e-^^^'' [1 + j^{t - u)] du 



infty 



Replace H[juJi) = Hie^°'^ 



H{joJi) + j / h{u)e-^'^^''ip{t - u) du 



Vo{t) = 3f?|e^('^i*+^) 



^ie^"i + 3 / h{u)e-^^^''ip{t - u) du 



1 + — — I h(u) je-^'('^i"+"iV(* - u) du 

Observe that je~-'('^i"+"i) = sin(a;iti + qi) + j cos(wiu + ai). Define 

hc{t) = h{t) cos(a;iii + ai) 
hs{t) = h(t) sin(cjiu + ai) 

Thus 



1 + 



hs{u)ip{t — u) du + j 



Hi 



hc{u)ip{t — u) du 
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Use 1 + a + j6 ~ (1 + a)e-'^ for small a and b 



o{t) = m e^^'^i'+^+^i^/Zie^'^iFi 



1 + 

Taking ai = 0, the equivalent filters are 

hs{t) 



{hs * ip){ty 




' .{he* ^p){ty 




Hi \ 


exp 


Hi \ 


} 



AM: 



H^ 



PM: 



H, 



From cos(a;) = ^ [e-^^ + e it follows that 



hcit) = - 



-jx 



h{t) 



Consequently 



2Hi 



H{j{io - LJi)) + H{j{io + uJi)) 



which is equivalent to Eq. H3.23() . ▲ 

Remark 3.5 Eq. ()3.23p states that the poles (and the zeros, if any) of 
H{s) appear in 7i{s) translated by ztwi. H{s) has (at least) a pair of 
complex conjugate poles on the left-hand half-plane at about iwo. A copy 
of them is present in TCis) translated close to the real axis. Another copy 
translated in the opposite direction is present in 'H(s), about ib2u;i. This 
second copy should be discarded because it results from the approximations 
taken. Anyway, questioning about the phase noise at a "close-in" frequency 
twice the carrier frequency is a nonsense. A 

Remark 3.6 If tJi ^ ujq, the resonator natural frequency, there results 
a zero on the real axis, in the negative half plane. This is still to be 
demonstrated with the formalism of this Section. ▲ 



Property The resonator converts phase noise into amplitude noise. The 
transfer function of this process is 



1 



2\H{jL0i)\ 



H{j{uj-uJi)-H{j{uj + uJi] 



(3.24) 
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Figure 3.9: Derivation of the oscillator phase-noise model (B) from the 
voltage-noise oscillator model (A). 



Proof. Proving Eq. (|3.23|) . we come across the AM filter ^ *. This 



Hi 

is the amplitude response to a phase modulation. Using sin(a::) 
j_ 
2i 



^ [e^^ — e ^^\, we get 



Kit) = 1 



h{t) 



Eq. (|3.24j) follows immediately. 



3.4 Phase noise in the oscillator loop 

Whereas phase noise can be analyzed with the generic model of Fig. 13.11 we 
turn our attention to the model of Fig. I3.9I B. specific to phase noise. The 
physical quantities of Fig. I3.9I B are the Laplace transform of the oscillator 
phase fluctuations. Thus the input signal ^{s) models the amplifier phase 
noise. Additionally ^{s) can be used to introduce the resonator fluctu- 
ations, or the phase noise of an external signal to which the oscillator is 
locked by injection. The amplifier gain is 1 because the amplifier repeats 
the input phase to the output. The filter transfer function /3(s) is replaced 
with -B(s), which is the Laplace transform of the resonator phase response. 
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Figure 3.10: Oscillator phase-noise transfer function TC(s). 



The phase-noise transfer function is 



^(s) 

n{s) = def. of n{s) 



n{s) 



1 



l-B{s) 



oscillator . 



(3.25) 
(3.26) 



Property The phase transfer function of the oscillator tuned at the exact 
resonant frequency is 



n{s) 



1 + sr 

ST 

t'^lo^ + 1 



exact natural 
frequency ujq 



(3.27) 
(3.28) 



Proof. Replace B{s) given by Eq. in Using = 

the result (|3.27|) follows immediately. The proof for \1-L{juj)\'^ is omitted. ▲ 



Leeson formula The oscillator phase noise spectrum 

2 



SM) 



P {2Qj 



(3.29) 
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where S^{f) is the spectrum of the phase fluctuation ip{t). 
Proof. From the definition of n{s) [Eq. ^71^ ]. it follows that 

S^H = |H(ja;)|2 5^M . 

We replace Eq. (|3.28j) in the above 

1 



Then, we replace ui = 27r/ and r = 22 = [Eq. (fXTIjl ] 



SM) 



1 + 



1 + 



1 



47r2/2 Q2 
P \2QJ 



Property The phase transfer function of the detuned oscillator, oscillating 
at the frequency loi = ujq + is 



2^2 



nis) 



detuned 
by O 



S{ST+ 1) 

(s + ^ - jn) (s + ^ + j^) alternate 
s (s + i) form 

T^oj^ + 2(t^ - n^T^)u;^ + (OV^ + 2^7V2 + 1) 



(3.30) 
(3.31) 

(3.32) 
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Proof. Replace B{s) given by Eq. I^TFIii in (HT^ . Use = 

^^^^ r{s+'--jn){s+'-+ jn) - (. + i + n^r) 
_ {sT + 1 - jnT){sT + 1 + jnr) 

~ {sT + l- jnT){sT + 1 + jnr) - (sr + 1 + 0^2) 

~ (sr + 1)2 + 172t2 - {sT + 1) - 02^2 

_ (sr + 1)2 + ^2^2 
~ s(sr+l) 

The proof for \7i{juj)\'^ is omitted. A 

Remark 3.7 Eq. (|i-{.H2j) extends the Leeson formula to the case of the 
detuned oscillator. A 

Figure Vd.lUl shows 'H(s) on the Laplace plane. When the oscillator 
is pulled out of luq, the real zero at s = —ujl splits into a pair of complex 
conjugate zeros at s = —ujLzizjQ, leaving a real pole at s = —ujl in between. 
The pole at s = in Fig. K-i.lOl is an ideal integrator in the time domain. 
This is the Leeson effect, which makes the phase diverge. At uj ^> ujl, 
7i{s) appears as a small cluster with equal number of poles and zeroes that 
null one another, for 7i{juj) is constant. The resonator is a flywheel that 
blocks the phase fluctuations, for the oscillator is open loop. Accordingly, 
the amplifier phase noise ^{s) is repeated as the output. 

The plot of |W(ju;)|2 (Fig.EHU) reveals that the noise response is a 
function of 0, and increases as O increases. This is best shown in terms of 
the normalized noise transfer function |'H„(ju;)|2, defined as 

1-^ r M2 def.of|?^„(ju;)|2 

\rin\J^)\ = -r ?5 — (3.33) 

|wo-)lLo Fis.Eia. 

The denominator \'H[juj)\^^Q is evaluated at = 0. |'H„(ju;)|2 has an 
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Oniega= 



Phase transfer function |H(jf)|»2 



tau=1 -> omega_L=1 
Omega = 

[0, 0.5, 1, 1.5, 2] omega_L 




0.2 0.4 



frequency f 



Figure 3.11: Phase- noise transfer function \7{{jLo)\'^ evaluated for r = 1. 



horizontal asymptote at |'H„(ja;)p = 1, for can not be integrated to infinity. 
After translating the horizontal asymptote to zero, we define 



M{n) = / [|W„(ja;)P - 1] dco 



(3.34) 
(3.35) 



A^(f2) is a parameter that describes the noise increase due to the detuning 
Q, which converges with actual resonators. Interestingly, it holds that 



vr 



M{^)='- n^T^ . 
4 



(3.36) 
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3.4.1 Pulling the oscillator frequency 

The oscillator can be pulled to a desired frequency in an interval around 
Wo, as shown in Section [1.2.11 

If the oscillator is pulled by inserting a static phase in the loop, the 
oscillation frequency is not equal to the resonator natural frequency. In this 
conditions it holds that ^ 0. Equation indicate that the oscillator 

phase noise increases. On the other hand, if ujq is changed by making the 
resonator interact with an external reactance, there is no reason for to 
be changed. Thus, the phase noise does not increase. The following cases 
are of interest. 

Quartz oscillator The quartz oscillator usually pulled by means of the 
external reactance, as in Figure fTTzl for Q is not affected. Yet, the additional 
loss introduced by the external reactance may lower the merit factor, and 
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in turn increase the resonator noise bandwidth. 



Low-noise microwave oscillators In low-noise schemes, hke the Galani 
oscihator (GBW"*" 84j . it often happens that the resonator can not be tuned, 
for a phase shifter is used to pull the oscillation frequency. There follows 
that 7^ 0, and additional noise is taken in. 

Pound-stabilized oscillator In the Pound scheme |Pou46j . the mod- 
ulation and detection mechanism ensures that the resonator is tuned at 
the exact natural frequency ujq by equating amplitude of the reflected side- 
bands. In this oscillator it holds that = 0. 



3.5 Spectrum of frequency fluctuation and AUan 



variance 



If the phase noise spectrum inside the loop is S^{f) = bo + 6-1/ ^ + . . ., 
the oscillator fractional frequency spectrum is 



^0 j;2 : J. , _bo_ I 
4g2 + 4Q2 / 



and the Allan variance is 
^2 



1/r 
terms 



+ l^bo- + 21n(2) ^ + . . . 



2 4g2 " r ' ' 4Q2 

The proof extends this result to higher order terms, 6-2/"^, etc. 
Proof (Spectrum) 



Syif) = ^SM) 



1 + 



P 4Q2 



derivative 
Leeson 



(3.37) 



(3.38) 
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+ 



^0 .2 \b-i bp b^i 1 



Proof (Allan variance) Match Eg. 13. 371 to Sy{f) = J2i b-iP, and iden- 
tify the terms 

Use Table O (page \^ . A 

Remark 3.8 For the reasons detailed in Section 11.41 terms steeper than 
6_2/~^ can not be present in the amplifier noise. They can be added to for 
completeness, because ^jJ(t) models all the phase fluctuations present in the 
loop. A 

Example 3.1 Calculate the Allan variance of a microwave DRO in which 
the resonator merit factor is Q = 2500, and the amplifier noise is S^{f) = 
10-^^ + 10-^7/ (white -150 dBrad^Hz, which results from / = 4 dB and 
Pp = -20 dBm, and flicker -150 dBradVHz at 1 Hz). Using Eq. iTTHHl) 

, 2x10-23 ^ ^ 4.47x10-12 
(j2(r) = + 5.55X 10-1^ cTy(r) = ^ + 7.45x10"^° 



▲ 



Chapter 4 



Noise in the delay-line oscillator 



The delay-line oscillator (Fig. I4.1I A) is an oscillator in which the feedback 
path is a mere delay r, for if the gain is ^ = 1 oscillation can take place at 
any frequency u;„ multiple of ^ in a range. A selector filter selects a given 
mode ujm = ^nT- by introducing attenuation at all other frequencies but 
LUm- The delay-line oscillator has the relevant property that the oscillation 
frequency can be switched between modes by tuning the selector, like a 
synthesizer. Of course, the pitch is Auj = 

We analyze the phase noise mechanism in the delay-line oscillator fol- 
lowing the same approach of previous chapter, based on the elementary 
theory of linear feedback systems. 



4.1 Basic delay-line oscillator 

The delay-line oscillator can be represented as the feedback model of Fig. 
I4.1I A. This is similar to Fig. l3.1l but for the resonator replaced with a delay 
line of delay r. The signal Vi{s) allows initial conditions and noise to be 
introduced in the system, as well as the driving signal if the oscillator is 
injection locked. We first study the voltage-noise transfer function, defined 
as H{s) = [Eq- (jS3)]- our system, it holds that H{s) = i_23(s) 
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A - Voltage space 



Vi(s) _ 



Vo(s) 







noise 








delay 









out 

v,;(s) 



oscill. 
out 



B - Phase space 



phase 
noise 



<I.(s) 



delay 
B(s) 



Figure 4.1: Voltage-noise (A) and phase-noise (B) oscillator model. 



[Eq. (|ll2j) ]. In the Laplace transforms, a delay r maps into e ^'^ . Thus 

^(^) = i-le-^r (Fig-OA) . (4.1) 

H{s) has an infinite array of poles that are solutions of "^(5) = 0, with 

V{s) = 1 - Ae-'^ = 1 - Ae-^e-'^^^ . 

The poles are located at s = s^, 

1 27r 

Sn = — In(^) +j — n integer n . (4.2) 
r r 

Figure shows the poles of H{s) on the complex plane. The poles are 

in the left half plane for A < 1, move rightwards as A increases, and reach 
the imaginary axis for A = 1. The derivative ^f- = indicates that a gain 
change 6 A cause the poles to move by 6sn = - Thus, the poles move 
horizontally by 6an = real, and vertically by jSu>n = ^ ^ if SA 

is imaginary. Yet for notation consistency, we will keep A a real constant, 
and introduce a separate function pe^^ (Sec. I4.2ll . 

A pair of imaginary conjugate poles represents a loss- free resonator that 
sustains a pure sinusoidal oscillation in the time domain if excited by an 
appropriate initial condition. Consequently, the delay-line oscillator can 
oscillate at any frequency lj„ = ^n. If more than one pole pair is excited, 
the oscillation is a linear superposition of sinusoids at the frequencies ujn- 
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Figure 4.2: Noise transfer function H{s), and phase noise transfer function 
fiis) on the Laplace plane. 



This is the Fourier series expansion of an arbitrary periodic waveform. The 
property of completeness of the Fourier series tells us that the delay-line 
oscillator can sustain any periodic waveform. 

Let us now choose one oscillation frequency, u>m- This means that only 
one pair of poles of H{s) is excited by appropriate initial conditions, at 
the frequencies uj±m = it^m. Having defined the carrier frequency oJm, 
we are now able to analyze the phase-noise transfer function, defined as 
TC{s) = ^7jy [Eq. (121251)] • The phase noise model of the oscillator is shown 
in Fig. I4.1I B. Once again, the amplifier gain is 1 because the amplifier just 
"copies" the input phase to the output. Thus it holds that 
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H(if)r2 




frequency f 



Figure 4.3: Phase noise transfer function evaluated for r = 1. 



which is the same as Eq. 1)3. 26|) . Yet, the phase transfer function is now 
B(s) = e~^'^ , independent of the carrier frequency. A proof can be obtained 
fohowing the approach of Section 13.21 Nonetheless, we observe that the 
delay line has infinite bandwidth, for the group delay is the same at all 
frequencies, and equal to the line delay r. Hence the phase, as well as any 
information-related parameter of the signal, take a time r to propagate 
through the delay line. Consequently it holds that 



By substituting B{s) = e 

n{s) 



B{s) = e-'^ . 
in 7i{s), we get 



1 



delay-line 
oscillator 



Tl{s) has an infinite array of poles that are solutions of 1 
Introducing the mode offset 



(4.3) 



(4.4) 



0. 



n 



m 



mode offset 



(4.5) 
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i.e., the integer frequency offset from the carrier, the poles (Fig. I4.21 B) are 
located at s = s^, 

= j— /i integer /i . (4.6) 

In the spectral domain, the square modulus of the phase noise transfer 
function is 

Wju^t = irj:^^ T • (Fig. USD (4.7) 

2 (1 — cos (jrj 

Interestingly, it holds that |W(i/)|^in = I- 



4.2 Mode selection 

A pure delay-line oscillator, with an infinite array of poles on the imaginary 
axis, is a nice exercise of mathematics with scarce practical usefulness. A 
mode selector is necessary, which ensures that that the Barkhausen con- 
dition A(3 = 1 is met at the privileged m only, thus at the frequencies 
oo±m. = ±^m, and that all the other poles Sn, n ^ ibm, are on the left- 
hand half plane (5?{s„} < 1). Thus we replace 

/3(s)^/3d(s)/J/(5) (4.8) 

in Eq. JUTJ), with 

/3^(s) = e-'^ delay line (4.9) 

/3^(s) = p(s)e-J'^(") selector filter (4.10) 

The function Pf{s), as any network function, has the following properties 

Pfis) = PKs*) i S^^!'?l = ^l^^i'?L T-"'' (4-11) 

^ ■^^■^ -1^ ' \ 9{/3/(s)} = -9{/3/(s*)} oddma; ^ ' 



p{s) = p{s*) the modulus p{s) is even in to (4-12) 

6{s) = —9{s*) the argument 9{s) is odd in lo (4-13) 
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These conditions are necessary for Pf{s) to be a real-coefficient function, 
or by extension a function whose series expansion has real coefficients of s, 
and ultimately for the inverse transform (i.e., the time-domain response) 
to be a real function of time. 

Introducing the selector ffiter, the transfer function [Eq. (|4.H) ] becomes 

= l-ApislMs)e-sr ■ (4-1^) 
The poles of H{s) are the solutions of 'D{s) = 0, with 

V{s) = 1 - Ap{s)e^'^^'^e-'^ (4.15) 
Splitting T>{s) into real and imaginary parts, we get 
Ape-"^ = 1 



and finally 



a = - \n{Ap) 

T 

e ^ 27r 

io = — mod — 

r r 



poles of i7(s) . (4.16) 



As a consequence, the modulus p affects only the real part a of the pole, 
thus the damping; the argument 6 affects only the imaginary part, thus the 
frequency. The poles of H{s) are located aX s = Sn 

Sn = -ln(^p) + j— n + j- poles of i?(s) (4.17) 
r r r 

The above Eq. H4.17|l will be used to analyze some relevant types of selector 
mter. 
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1 1 



Figure 4.4: Amplitude-only selector filter. 

4.2.1 Amplitude-only filter 

This selector filter is defined as 

^ UJ ^ ±LOm 

, N I (n = ibm) , __ 

(3f{s) = { \ > (Fig. nan. (4.18) 

eisewnere 

[n 7^ ±m) 

Such filter attenuates by a factor e ^ 1 all the signals but those in the 
vicinity of icjm. Needless to say, this filter is an abstraction, which can 
not be implemented with real-world devices. 

The poles of H{s) are the solution of T){s) = 0, that is, 

(n = ±m) ^ ' 

Hence the poles (Fig. I4.5I A) are located at s = s„, 

Mn(^)+i^n n = ±m (w = itJ^) 
\ In(eA) -h j^n n / ±m 



(4.21) 



4.2.2 Phase-only filter 

The phase-only filter is defined by 

/3^(s) = e^''^snM^ (4.22) 
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Figure 4.5: Transfer function H{s) of the delay-line oscillator with a selector 
filter /3/(s) inserted in the oscillator loop. 



This filter shifts the signals by 9 for cu > 0, and hy —9 for u; < 0. The 
signum function sgn(a;) is necessary for the condition f3f{s) = f3*f{s*) to be 
satisfied. Once again, this filter is an abstraction. 
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By replacing the above /3/(s) in H{s), we get 

1 



H{s) 



1 — Ae^ sgn(a;) 

The poles are the solution of D = 0, where 

V{s) = 1 - Ae^'«§°H^e-^^ 
By splitting the modulus and the argument of D, we find 
Ae-"^ = 1 

sgn(a;) 9 — ujt = mod 2tt . 
There follows that the poles (Fig. I4.5I B) are located at s = Sn 
1 



(4.23) 



■ HA) + J 



2tt . . 
— n H — sgn(nj 
r r 



integer n 



(4.24) 



In conclusion, the real part of the poles is not affected, for the gain condition 
is still A = 1. On the other hand, the phase 9 results in a frequency shift 



Ao; = — for w > 
r 



and 



for w > . 



(4.25) 



4.3 The use of a resonator as the selector filter 

The practical reason to choose a resonator as the selector filter is that it 
is simpler and easier to tune, as compared to other types of filter. Let us 
assume that the resonator is tuned at the exact frequency ujm = of 
our interest, among the permitted frequencies of the delay-line oscillator. 
Thus, the feedback function is 



Pis) = Pd{s)l3fis) 



[Eq. ^] , (4.26) 
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with 

^'^^^^-'^ [Eq. gSI)] ^ ' 

a ( \ ^rn s filter , . 

^f^'^ - -Qs^ + ^s + ul [Eq. ^^-'^^ 

It is important to understand the role of delay line and of the filter. As a 
consequence of the Leeson effect, stable oscillation frequency requires high 
slope ^ arg/3(ja;), i.e., long group delay. Nonetheless, the oscillator tracks 
the natural frequency of the feedback network, which is sensitive to the 
environment parameters, chiefly the ambient temperature. When (3d and 
(3f are cascaded, the frequency fluctuations are weighted by the phase slope 
^ arg /3{ju>) because the Barkhausen condition requires arg P{juj) = 0. The 
appropriate choice for the delay-line oscillator is therefore 



dco 



(4.29) 



which means that the oscillator frequency is stabilized to the delay line, 
while the resonator serves only as the mode selector. The main technical 
reason to choose a delay-line oscillator is the availability of a delay line that 
exhibit much higher stability than the resonator. Accordingly, the resonator 
merit factor must be sufficiently low to prevent the resonator from impairing 
the oscillator stability. Another reason for the condition (|4.29|) is that it 
makes easy to switch between the delay-line modes, allowing the resonator 
to be tuned in a narrow range around the exact frequency. 

By replacing /3(s) = (3d{s)[3f{s) in H{s) = iz^xpTs)'' under the assump- 
tion that the amplifier gain is A = 1, we find the oscillator transfer function 

This function has a pair of complex conjugate zeros at s = and s = s*^ 
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and a series of complex conjugate poles. As a consequence of the condition 
H4.29() . the resonator bandwidth ^ is large as compared to the mode pitch 
This means that in a frequency range J- around Wm, and also around 
—iOm, the resonator dissonance x is such that Qx ^ 1- The range 
contains a few modes a;„ around Wm- In ^ the resonator function /3j(ju;) 
is about constant and close to 1. Thee follows that in 

1. The poles of H{s) are chiefly determined by the oscillation of Pd{s)', 
they are close to j^n, which is the result already obtained without 
the selector filter. 

2. An approximate solution for the poles is found by replacing /3/(s) 
with Pfijuj) = p(cj)eJ'^('^) 



X 



— arctan (Qx) 



UJ 

UJ. 



m 



UJ 



modulus 

phase 
dissonance 



3. From Qv < 1 it follows that y ~ 2"^^^ 
6{uj) can be further approximated as 



Consequently p{uj) and 



p{uj) 
e{uj) 



1-2Q 
-2Q 



U! — UJr, 



UJ. 



m 



Thus, the poles of H{s) are found by inserting the above p{uj)\i^=uj„ and 
9{^)\ui=Lu„ in Eq. (|4.17jl . with A = 1. The poles are at s = s„ = (7„ + jujn 



1 



■In 



T 

27T 



UJr. 



UJn 



2Q2 



UJn 



UJ.r, 



UJn 



e 2-K 



2Q UJr. 



(4.31) 
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hence 



2Q'^ f IVn - UJm\'^ , .2vr 2QuJn-UJm 

+ j — n-j , (4.33) 



Additionally, the pols can be given in terms of the mode offset ^ = n — m 
[Eq. (|4.5)) ] through the simple property that ^"J"*^'" = Hence 

,„ = -?e!(ii)V,^„-,^ii. (4.34) 

r \mJ T T m 

Figure shows the poles of H{s), located on a horizontal parabola cen- 
tered at ujm- The pole frequencies are shifted by — with respect to 
the exact mode of the delay line. This frequency shift results from 
the off-resonance phase of the filter. The negative-frequency part of the 
Laplace plane follows by symmetry. 

Interestingly, every pair of complex conjugate poles of H(s) can be 
seen as a separate resonator [Fig. lA.ll and Eq. ()A.1|) ]. whose merit factor is 
enhanced by the positive feedback. The pole pair selected by the resonator, 
denoted with the subscript m, is onto the imaginary axis, for the equivalent 
merit factor is Qm = oo. The other pole pairs, denoted with the subscript 
n, yield a merit factor 

Qn-'^(-Y (4.35) 

or equivalently 

Qn 



4Q2 V/i 



vrn nn 



2Q2 V /i . 

In the presence of noise, these poles are excited by random signal. Thus, 
the radio-frequency spectrum shows a series of sharp lines at frequency lo^ 
off the carrier frequency tOm, on both sides because /i takes positive and 
negative values. These lines are easily mistaken for competitor oscillation 
modes, which they can not be^. The relevant differences are that random 

^Mode competition and multimode oscillation do exist, yet under different hypotheses. 
This behavior is commonly found in lasers. 
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Figure 4.6: Transfer function H{s) of the delay-line oscillator with a res- 
onator as the selector filter. 



excitation causes incoherent oscillation, and that such oscillation is damped. 
With the /i = pair, coherent and incoherent oscillation coexist. This is 
the Leeson effect. 

With the operating parameters of actual or conceivable oscillators, 
the equivalent merit factor Qn is so high that the noise transfer function 
\H{jLi!)\'^ [Ea. l4.3Uj can hardly be plotted. Examples are provided in Section 

EM 



92 



E. Rubiola The Leeson effect 



February 2, 2008 



Remark 4.1 (Multi-pole filters) In the domain of telecommunications, 
multi-pole passband filters are widely used, which exhibit nearly flat re- 
sponse in a bandwidth. The experimentalist may be inclined to use such 
filters as the mode selector because of their easy commercial availability, 
and out of the common belief that the flat response is the most desirable 
one for a bandpass filter. But the maximally flat filter is the worst choice 
for an oscillator. In fact, with a flat filters the 7^ poles are maximally 
close to the imaginary axis. There follow that 

1. selecting the desired oscillation frequency is maximally difficult, 

2. mode jump can occur, 

3. close-in noise peaks are maximally high. 



4.4 Phase response 

Let the delay-line oscillator oscillate at the frequency ojm = thanks to 
a resonator in the feedback loop tuned at the exact frequency tOm- Under 
this condition, we analyze the phase-noise transfer function 'H(s) = 



Introducing the resonator as the selector filter, the feedback function i 



A 



[Eq. dS^ni)], i-e. 




(4.36) 




(4.37) 



with 




delay 



(4.38) 




l + STf 



filter 



(4.39) 
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Figure 4.7: Phase-noise transfer function 'H(s) of the delay-hne oscillator 
with a resonator as the selector filter. 



The subscript 'd' and '/' are introduced in order to prevent confusion be- 
tween the delay of the line (r^) and the relaxation time rj = ^ [Eq. ()A.14|) ] 
of the resonator. 

Expanding the phase-noise transfer function, we find 

n{s) = ^-±111 (4.40) 

Tl{s) has a real zero at s = = — f^, a pole at s = 0, and a series of 
complex conjugate poles in the left half-plane. The location of these poles is 
best found by searching for the zeros of 1—B{s), without using the simplified 
form (|4.4()j) . As a consequence of the design choice Tf = the poles 
are expected to be close to the position already found in the absence of the 
selector, that is, jf"/^- In other words, in a region around the origin that 
contains a few pole pairs, the selector has a small effect on the modulus 
and on the phase of B(s). This is exactly the same situation already found 
searching for the poles of H(s), transposed from ±uJm to the origin. 

In order to the find poles, solution of 1 — B{s) = 0, we first focus on 
the low-pass phase filter Bf{s). The poles are expected to be close to the 
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imaginary axis. Thus we can approximate Bf(s) ~ Bf{jio), and write in 
polar form 

Bf{s) = p(cj)e^'''(^) 
p{u;) = 



1 



1 + uj^t] 



9{oj) = — arctan (ujTf) . 

Then we replace the low-pass time constant Tf as tj = and we expand 
for iOTf <^ 1 



e{uj) = -2Q — 
The equation 1 — B{s) = can be rewritten as 

and further split into modulus and phase. The poles are the solutions of 

p(a;)e~'^^'* = 1 modulus 
9{uj) — ujT(i = mod 27r phase . 



The modulus condition yields 

OJ, 



- 1^ 



aTd = In 1 - 2Q 



m 

2 



12 ' 
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The phase condition yields 

U! 

LOTfi + 2Q = mod 2tt 

2tt 2Q u: 

u! = — fi integer /i . (4.42) 

The poles are located at s = = o"^ + ju^, which results from Eq. (|4.4H) 
and (jO^ 

.^._?5!4+,?^,_,?0^ (4,43) 

Td Td Td UJm 

Additionally, by replacing ^ = /i we find 

2Q2.^x2 2^ 2Q ^i 
Sf, = — +j— /i . 4.44) 

Td \mJ Td Td m 

Figure 14.71 shows the phase- noise transfer function 'H{s) on the Laplace 
plane. The pole in the origin represents a pure integrator in the time 
domain, which causes the Leeson effect. The other poles are on a horizontal 
parabola centered in the origin, at a small negative distance a oc — from 
the imaginary axis. It is easily seen by comparing Fig. 14.71 to Fig. lA.ll that 
each complex conjugate pair is equivalent to a resonator of merit factor 

which can be rewritten as 

because Un = —a. 

From Eq. (|4.4()|) . after some tedious algebra it is proved that 

\n{3^)? = ^^ -, r (Fig.lIHI). (4.47) 

2 — 1 cos[u;Td) + ^ T^ + 2ujTf sm[LOTd) 
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- Phase transfer function |H(f)|»2 

1e+07. 

" parameters 



nu_m = 10 GHz 
- selector Q = 2000 




o 

frequency, Hz 

.1e-1-| 1 1 — I I I I I 1 1 1 1 — I I I I 1 1 1 1 1 — I I I I 1 1 

1 .1e4 .1e5 1e+05 1e+06 
O 



Figure 4.8: Phase transfer function \'H{jf)\'^ evaluated for = 20 /is, 
m = 2x10^ (thus Um = 10 GHz), and Q = 2000 (data of Example EH) . 



\7i{ju))\'^ shows a series of sharp peaks at u; = or / = ^fi, integer 

fi, around which phase noise taken in is enhanced. The peaks derive from 
the poles of 'H{s)'^ close to the imaginary axis, but not on the axis, for 
they are not the signature of competing oscillation modes. The minima of 
\7i{ju>)\'^ are located at u; ~ or / ~ // = ^, |, |, . . .. Interestingly, 
at these minima it holds that \7i{jui)\'^ ~ ^. Finally, at low frequencies 
(w < 1^ < 1^) |W(jw)p is approximated as 



4.5 Close-in noise spectra and Allan variance 



The phase noise spectrum is 

sM) = \njf)\'s^{f) 



(4.49) 
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At low frequencies, it holds that 

SM) j-2 SAf) '-«^ (4-50) 

^S^if) nseu^ = ^ (4.51) 



^ ^m^o J_ , J_ amplifier 

47^2^2 j2 ^ ^^2^2 p noise ^ > 

The spectrum of the fractional frequency fluctuation y is found using 

SyU) = ^S^U) 

Sy{!) = ^,S,i!) (4.53) 

By replacing Sjp{f) = bo + 6-i/~^ and matching the above to Sy{f) = 
Y^ihiP, we find 



Using Table imi the Allan variance is 



+ 7^^ + r^21n(2) 



4.6 Examples 



Example 4.1 Photonic delay-line oscillator We analyze a photonic 
delay-line oscillator based on the following parameters 

Td = 20 ^s 4-km optical fiber, refraction index of 1.5 

l/^d = 50 kHz mode spacing 

Urn = 10 GHz selected frequency 

m = 2x10^ mode order, equal to UmTd 

Q = 2x10^ selector filter (tunable microwave cavity) 

6o = 3.2x10"^^ white phase noise, —135 rad^/Hz 

6_i = 10"^° flicker phase noise, -100 dBrad^Hz at / = 1 Hz. 
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Figure EHl shows the phase noise transfer function \TL{jf)\'^. A series of 
spectral lines due to noise is evident at / = = 50 kHz and multiples. 
The following table shows the relevant parameters of the resonance due to 
the neighboring modes of the delay line. 







Qx 




Qn 
















OO 


OO 


±1 


5x10^ 


0.02 


-10 


3.14x10^ 


15.7x10^ 


±2 


10x10^ 


0.04 


-40 


785x10*^ 


7.85x10^ 


±3 


15x10^ 


0.06 


-90 


349 xlO^ 


5.24x10^ 


±4 


20x10^ 


0.08 


-160 


196x10^ 


3.93x10^ 


±5 


25x10^ 


0.10 


-250 


126x10*^ 


3.14x10^ 



The reader should observe the high equivalent merit factor of the phase- 
noise response, which explains the sharp resonances of Fig. 14.81 

The amplifier corner frequency is fc = 3.16 kHz, as results from 
6_i/~^ = Bq. The estimated phase noise spectrum is shown in Fig. 14.91 
Restricting our attention to the white and flicker frequency noise, the os- 
cillator phase noise is 

2><l(r^ 6.3x10-3 
^ifU ) — -p I js • 

The spectrum of the fractional frequency fluctuation is 

and the Allan variance 

-|^Q-26 IQ-'^'' 
al{T) = + 8.8x10-23 thus cjy(r) = — — + 9.4x10-^^ 



Remark 4.2 The numerical values of Example 14 . 1 1 are inspired to the ref- 
erences |E1V1()3| and |YlV196j . ▲ 
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Figure 4.9: Expected phase noise spectrum S^{f) for a photonic delay-hne 
oscillator. Parameters: r = 20 //s, m = 2x10^ (thus Um = 10 GHz), and 
Q = 2000. Data refer to Example O 



Example 4.2 SAW delay-line oscillator We analyze a SAW delay-line 
oscillator based on the following parameters 

15 mm SAW, sound speed 3 km/s 
mode spacing 
selected frequency 
mode order, equal to I'mT 
selector filter {LC filter) 
white phase noise, -145 radVHz 
flicker phase noise, —100 dBrad^/Hz at / = 1 Hz. 

The amplifier corner frequency is fc = 31.6 Hz, as results from b-if~^ = 
bo. Figure EHl shows the expected phase noise S^{f). The neighboring 
resonances show up at / = ^ = 200 kHz and multiples. The following table 
shows the relevant resonance parameters. 



T = 5 ^s 
l/r = 200 kHz 
= 900 MHz 
m = 4500 
Q = 80 

bo = 3.2x10-1^ 
6„i = 10-13 



100 



E. Ruhiola The Leeson effect 



February 2, 2008 




Figure 4.10: Expected phase noise spectrum S'(^(/) for a SAW delay-line 
oscillator. Parameters: r = 20 /is, m = 2x10^ (thus Vrn = 10 GHz), and 
Q = 2000. Data refer to Example Ol 





W-i^m\ 


Qx 




Qn 
















OO 


OO 


±1 


2x10^ 


0.035 


-126 


2.24x10^ 


4970 


±2 


4x10^ 


0.071 


-506 


5.59x10^ 


2485 


±3 


8x10^ 


0.107 


-1138 


2.49x10*^ 


1657 


±4 


10x10^ 


0.142 


-2023 


1.40x10^ 


1243 


±5 


12x10^ 


0.178 


-3160 


8.95x10^ 


994 



The estimated phase noise spectrum is shown in Fig. 14.101 Focusing 
our attention to the white and flicker frequency noise, the oscillator phase 
noise is 
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The spectrum of the fractional frequency fluctuation is 

5,(/) = 3.96x10-2^ + 1:^^^^^, 

and the Allan variance 

2x10-24 ^ _ _ _ 1.4x10-13 

7^ 



cj2(r) = + 1.7x 10-22 thus cJj^(t) = y= + 1.3x10-^2 



Remark 4.3 The numerical values of Example 14.21 are inspired to the eu- 
ropean GSM standard for mobile phones, which operates in the 900 MHz 
band with a channel spacing of 200 kHz. By tuning the selector filter, the 
SAW oscillator can be used as a frequency synthesizer that works at the 
appropriate frequency, with the appropriate step. A 



4.7 Phase noise in lasers 

Figure B.lll sketches the optical system and the equivalent electric circuit 
of a laser. T and R are the mirror transmission and reflection coefficients, 
A is the round-trip gain, and r is the round-trip transit time. Gain and 
transmission coefficients refer to electric field, thus to voltage in the equiv- 
alent circuit of Fig. l4.ljl B. We have to prove that the laser is equivalent to 
the feedback scheme of Fig. 14.11 for all the framework of Sections 14.31 14.41 
and 14.51 applies. 

At the reference plane, the input signal (electric field) is added to the 
signal amplified and fed back after delay r. Without loss of generality, we 
drop the trivial factor Ti. This means that the true input signal is Vi on 
the reference plane, instead of the external signal Then, we take the 
output signal Vo at the reference plane, instead of after the mirror 2. 
The transfer function is defined as 

H{s) = [def. of His)] . (4.54) 
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A - Optical system 



B - equivalent circuit 



reference 
plane 



v; T 




reference 
plane 



v; 



V: 




v„ 



out 



A' 



(out) 



delay 
X 




mirror 1 mirror 2 

Figure 4.11: A laser can be regarded as a feedback system. 



By breaking the circuit at the reference plane and equating the left-hand 
signal to the right-hand signal, we find 

Vois) = Vi{s) + Vo{s) [1 - A'{s)/3d{s)RiR2)] 

where (3(i{s) = e"*"^** is the delay. Hence 

^^'^ ^ 1 - A'{s)f3d{s)RiR2 

We observe that in actual lasers the gain A' of the active medium is a 
smooth function of frequency, wide as compared to the mode spacing of 
the optical cavity, which shows a maximum centered at the laser frequency. 
This is the same analysis/design criterion that we have adopted in Section 
14., S( introducing the mode-selector filter. Thus 
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1. We move the frequency dependence from A'{s) to a new filter function 
(3f{s). The latter is normalized for the maximum to be equal one, at 
the laser frequency. 

2. We approximate wide and smooth around the laser frequency, 
with a resonator. 

3. We introduce the gain A that accounts for the maximum and 
for the two reflection coefficients Ri and i?2- 

After these manipulations, the laser transfer function takes the form 



' ' 1 - APf{s)(3dis) ' 

which is equivalent to Eq. (|4.3()|) . Needless to say, the Barkhausen condition 
still holds, for the loop gain must be equal to one as a result of some 
saturation phenomena. 

4.8 Saturation and multimode oscillation 

Multimode oscillation requires that the Barkhausen condition Af3 = 1 is 
met at more than one frequency. This is of course possible. 

Single-mode oscillation In the amplifiers of the first type, when the sig- 
nal exceeds a threshold the gain drops more or less smoothly, yet uniformly 
in the pass band (Fig. I4.12I A). This is always the case of the oscillators 
in which the power is controlled in closed loop by measuring the output 
power and by acting on a variable-gain element. This behavior is also 
typical of saturated electronic amplifiers, either vacuum or solid state. In 
this case, distortion reduces the gain and squeezes power into higher har- 
monics, at frequencies multiple of ujq. Additionally, this is found in the 
quantum amplifiers (masers and lasers) in which amplification is due to a 
single pumped level, or to a cluster of levels that act as a single one be- 
cause of the linewidth. Oscillation sinks atoms from the pumped level, for 
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A - single gain 



B - gain cluster 



arg(A|3) =0 
largest gain 

small signal \ 



arg(A|3) 



large signal 




lAPI 



large signal 

|A|3|=1 arg(A|3)=0 
small n " 
signal Y 




arg(A|3) 



CO 



Figure 4.12: Different types of saturation in amplifier. 



the active population is reduced. With these amplifiers, the gain condition 
can only be met at one frequency. Generally, the oscillation frequency is 
that at which the small-signal gain is the highest among those at which the 
loop phase is zero. In the example of Fig. I4.12I A. u>i is not in the mode 
competition because of the insufficient the small-signal gain, and wins. 

Multimode oscillation In the amplifier of the second type, found in 
some lasers, the bandwidth is due to a cluster of energy levels narrow enough 
for power to sink population selectively (Fig. l4.1^ B). When oscillation takes 
place at a given frequency, the active population decreases only in the 
neighbors. Thus oscillation can also rise at the other frequencies at which 
the phase condition is met and the small-signal gain is greater than one. 
In the example of Fig. I4.12I B. the small-signal gain is insufficient at wi; 
oscillation takes place at 102, 1^3, u;4, and ws. 



Appendix A 



Resonator model and parameters 



A.l Laplace plane 

Around the resonant frequency, real world resonators are well described 
by a linear second-order differential equation with constant coefficients. 
Accordingly The Laplace transform of the impulse response is a rational 
function of the form 

which can be rewritten as 

Hr{s) = ^ r + (A.2) 

with 



Wo 
2Q 



(A.3) 



^l^^l = ^l- (A.5) 
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Figure A.l: Resonator transfer function in the frequency domain and on 
the Laplace plane. 



Figure lA. II shows the resonator transfer function in the frequency domain 
and on the Laplace plane. 

A. 2 Frequency domain 

The frequency response is easily found by replacing s = joo in the Laplace 
transform. Defining the dissonance x ^-s 

X = — — ~ dissonance (A. 6) 

we find the following relevant relationships 

Hr{juj) = ^ _^ -Q frequency response (A.8) 

R{M = 1 , L 2 real part ^{Hr{ju;)} (A.9) 
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Qx 



1 + Q2^2 

1 

— arctan Qx 



imag. part ^{Hr{jijj)} (A. 10) 

modulus \Hr{juj)\ (A. 11) 

argument arg [ffr(ju))] (A- 12) 



A. 3 Time domain 



In the time domain, the impulse response of the resonator an exponentially 
decaying sinusoid of the form 

Vo{t) = cos {uQt) e"*/^ , (A. 13) 

where the relaxation time r is related to the other resonator parameters by 

1 relaxation ^^^^^^ 



relaxation 
time 
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